ii 


TABLE  OF  CONTENTS 


Page 


Acknowledgements  ..........  .  iv 

List  of  Tables  ........  .  v 


List  of  Figures 


Nomenclature  .  .  .......  . 

V 

CHAPTER  I  INTRODUCTION 

1.1.  Reasons  for  Studying  Spherical  Shockwaves.. 

1.2.  Physics  of  Cavity  Expansion  due  to  Explosive 

Impact, .  .  .......  . 

1.5*  General  Nature  of  Shock  Waves  .  .  . 

1.4.‘  Statement  of  Problem  and  Scope  of  Thesis:, 

CHAPTER  II  BRIEF  REVIEW  OF  PAST  WORK-  - 

2.1.  Elasti c -Plastic  Effects,  .  . 

2.2.  Similarity  Solutions  of  Equations  of  Gas  Dynamics. 
2-3-  Self -Similar  Solutions  of  the  Problem  of  Cratering 

Due  to  Hypervelocity  Impact,.  ....... 

2. 4T  "Numerical  Integration  of  the  Differential  Equation 
of  Spherical  Blast  Waves  in  a  Gas  and  Hydrodynamics 
of  Hypervelocity  Impact;  ......... 


CHAPTER  III  SELF -SIMILAR  SOLUTION  FOR  PROBLEMS  WITH  SPHERICAL 
SYMMETRY  -  - 


3.1- 

3.2." 


Theory, 

Solution  of  the  Governing  Differential  Equation 


of  Self -Similar  Motion,' 


CHAPTER  IV  APPLICATION  OF  THE  THEORY  OF  SELF -SIMILAR  MOTION 
TO  THE  PROBLEM  OF  EXPANSION  OF  A  SPHERICAL  CAVITY 
IN  A  METAL  DUE  TO  EXPLOSIVE  IMPACT  -  - 


4.1.  Description  of  Model  and  Assumptions, . 

4.2.  Equations  of  State  of  Metals^  . 

4.3.  Solution  of  the  Blast  Wave  Problem  in  Aluminum., 

4.4.  Blast  Waves  in  Other  Metals, . 

4.5 "Consideration  of  the  Non-Similar  Aspects  of  the 

Blast  Wave  Problem 

4.6.  Summary  and  Conclusions  .  . 

L  n  1 


CHAPTER  V 
5-1. 

5-2. 

5-3. 


EXPERIMENTAL  PROCEDURE  AND  RESULTS; 
Descriptions  of  Specimens  and  Experimental 

Arrangement  . 

Test  Results  . 

Discussion  and  Conclusion  . 


vi 

vii 

1 

2 

2 

4 

6 

7 

7 

8 


10 

21 


23 

23 

31 

44 

49 

54 


58. 

61 

63 


1 


ill 


Page 

CHAPTER  VI  CLOSURE 

6.1.  Summary  and  Conclusions  . .  JO 

6.2.  Suggestions  for  Further  Research .  73 

■  iJ-  •  \ 

APPENDIX  A  'SINGULAR  POINTS  OF  ORDINARY  DIFFERENTIAL  EQUATIONS  j  75 
APPENDIX  B  NUMERICAL  PROGRAM-;  C  79 

BIBLIOGRAPHY  '  85 


iv 


Acknowledgements 

The  work  of  this  paper  was  sponsored  by  the  Ballistic  Research 
Laboratories,  Aberdeen  Proving  Ground,  Their  support  is  gratefully 
acknowledged. 

Acknowledgement  is  also  due  to  0.  T-  Johnson  of  Ballistic  Re¬ 
search  Laboratories,  for  performing  the  experiments  and  A.  Hoffman 
end  R.  E.  Shear,  also  of  Ballistic  Research  laboratories,  for  many 
helpful ' suggestions . 


vi 


LIST  OF  FIGURES 

Figure  Page 

1  Hugoniots  of  a  Polytrqpic  and  Quasi -Polytrppic  30 

Medium 

2  Variation  of  the  Functions  f^,  f 2,  f ^  with  c/C  33 

3  Solution  Curve  for  Aluminum  35 

4  Hugoniots  of  Four  Metals  3 6 

5  R-T  Diagram  for  Shock  Region  for  Aluminum  Sphere  40 

6  Cavity  Expansion  and  Particle  Trajectories  42 

7  Cavity  Expansion  in  Four  Metals  45 

8  Velocity  Decay  on  Cavity  Surface  for  Four  Metals  46 

9  Pressure  Decay  on  Cavity  Surface  for  Four  MetalB  47 

10  Relative  Position  of  S4  and  Shock  Front  for  50 

Variable  7 

11  Smaller  Sphere -Pre  Shot  59 

12  Experimental  Arrangement  62 

13  Response  of  Capacitance  Gage  for  Smaller  Sphere  65 

14  Response  of  Capacitance  Gage  for  Larger  Sphere  67 

15  Smaller  Sphere -Post  Shot  68 

16  Behavior  of  Solution  Curves  for  Variable  7  69 


2 


1.2.  Physics  of  Cavity  Expansion  Due  to  Explosive  Impact 

Consider  a  detonation  Initiated  at  the  center  of  a  spherical 
charge  of  high  explosive  surrounded  by  a  ductile  metal.  A  wave  races 
through  the  explosive  and  imparts  a  sharp  impact  to  the  surrounding 
material.  A  shock  wave  is  formed  in  the  metal  raising  the  the  pressure 
to  the  order  of  magnitude  of  500  kilcbars.  Under  theee  extreme  pres¬ 
sures  the  metal  in  a  zone  adjacent  to  the  explosive  behaves  as  a  com¬ 
pressible,  inviscid  fluid.  The  shock  wave  propagates  outward  causing 
large  plastic  deformation,  but  decays  due  to  spherical  divergence,  that 
is  the  pressure  acts  over  an  increasing  area,  into  a  stress  wave. 

Meanwhile  in  the  gas  there  is  a  shock  wave  propagated  back  toward 
the  center  of  symmetry  which  is  reflected  again  causing  the  gas  to  pul¬ 
sate  until  the  pressure  is  sufficiently  reduced. 

Vibration  or  oscillation  also  takes  place  in  the  solid  and  for  a 
finite  body  the  problem  of  interaction  from  reflected  waves  arises 
Thus  it  becomes  evident  that  the  complications  produced  by  the 
extreme  pressures,  the  varying  mechanical  behavior  of  the  metal,  short 
time  of  the  phenomena  as  well  as  the  interaction  between  gas  and  solid 
pose  a  difficult  problem  in  the  study  of  shock  waves  in  solids,  if 
attempted  in  its  full  generality.  An  analytical  theory  can  be  formu¬ 
lated  and  carried  through  only  by  isolating  the  significant  factors  in¬ 
fluencing  the  behavior  of  the  metal.  This  will  be  done  in  later  chapt¬ 
ers  of  this  work. 

1.5.  General  Nature  of  Shock  Waves 


The  concept  of  a  shock  wave  is  often  misunderstood  even  among 
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CHAPTER  II 

BRIEF  REVIEW  OF  PAST  WORK 
2.1.  Elastic -Plastic  Effects 

The  most  significant  work  published  on  the  problem  of  expansion  of 

a  spherical  cavity  in  which  elastic  and  elastic -plastic  effects  are  em- 

( l)* 

phasized  was  that  by  Hopkins  .  This  article  is  a  survey  of  the  lit¬ 
erature  and  includes  comments  on  most  of  the  significant  work  associat¬ 
ed  with  the  cavity  expansion  problem  to  date.  Solutions  are  given  for 
the  elastic  problem  using  the  classical  approach  of  integral  transforms. 
It  is  of  Interest  to  mention  that  in  the  discussion  of  elastic -plastic 
waves,  Hopkins  emphasizes  that  due  to  the  inhcmogeneity  of  the  dis¬ 
placement  equation  of  equilibrium  the  plastic  disturbances  are  not 
propagated  at  a  definite  velocity.  Only  that  part  of  the  plastic  wave 
derivable  from  the  plastic  displacement  potential  travels  at  the  plas¬ 
tic  wave  velocity. 

It  is  also  shown  that  even  the  relatively  simple  case  of  elastic- 

plastic  effects  In  a  perfectly  plastic  or  Reuss  type  material  that  the 

location  of  the  elastic -plastic  boundary  is  not  obtainable  in  closed 

(2) 

form.  An  alternative  procedure  suggested  by  Hunter  is  to  specify 
the  location  of  the  interface  and  then  determine  the  corresponding  wave 
functions  and  applied  pressure.  The  results  obtained  from  this  ap¬ 
proach  were  in  most  cases  found  to  be  unrealistic  or  to  place  such 
severe  limitations  on  the  applied  pressure  as  to  be  of  little  use. 

♦Figures  in  parenthesis  appearing  as  superscripts  refer  to  references 
listed  in  Bibliography. 
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Davids  and  Huang  Ky  were  apparently  the  first  to  approach  the 
problem  of  cratering  using  the  similarity  or  progret  fng  wave  theory. 
They  applied  the  theory  of  self -similar  flow  using  a  perfect  gas  type 
equation  of  state  for  the  target  material  and  an  energy  criterion  for 
establishing  the  crater  depth.  They  assumed  spherical  symmetry  in  the 
problem. 

Eae  and  Kirchner  formulated  the  problem  more  generally  but 

succeeded  in  getting  closed  form  solutions  for  the  spherically  symmet¬ 
ric  case  only.  They  also  pointed  out  some  restrictions  on  the  equation 
of  state  and  suggested  ways  of  approaching  the  non-similar  aspects  of 
the  problem.  A  material  strength  criterion  was  used  to  predict  crater 
depth . 

2.4.  Numerical  Solution  of  the  Problems  of  Blast  Waves  In  a  Perfect 
Gas  and 


Numerical  solutions  of  the  problem  of  Btror.g  blast  waves  in  a  per¬ 
fect  gas  have  appeared  in  the  literature  since  about  1954.  For  the 
most  part  the  approaches  are  all  similar  but  there  are  seme  distinct 
differences . 

Goldstine  and  Von  Neumann  solve  the  problem  of  a  strong  point 

explosion  in  an  ideal  gas.  The  conservation  laws  of  mass  and  momentum 
are  combined  and  the  resulting  differential  equation  solved  by  finite 
differences.  They  employ  an  iterative  scheme  to  satisfy  the  Rankine- 
Hugonlot  conditions  and  use  this  as  a  check  on  the  accuracy  of  their 
solution. 


Brode  claims  to  be  following  the  method  of  one  of  the  previous 
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CHAPTER  III 

SELF- SIMILAR  SOLUTION  OF  PROBLEMS  WITH  SPHERICAL  SYMMETRY 
3*1*  Theory 

The  theory  of  self-similar  motion  was  originally  conceived  for  the 
solutions  of  problems  in  gas  dynamics.  However,  subsequent  investiga¬ 
tions  have  led  to  the  admission  of  a  broader  class  of  problems. 

(a)  Basic  Equations  of  Spherical  Flow 

The  equations  for  spherical  flow  of  a  medium  in  Eulerean  form 
expressing  the  conservation  of  momentum, mass  and  energy  are 


du 

5t 


+  u 


0 


(3.1) 


dp 

5t 


dp 

u  5F  + 


,du 

(3F  + 


2u 


(3-2) 


S+’5t<F>-° 


where 


D  d  d 

Dt  =  5t  +  u  3r 


(3-3) 


and  e  =  f  (p,p)  is  the  form  of  the  equation  of  state. 

The  Rankine-Hugoniot  equations  expressing  conservation  of 
mass,  momentum  and  energy  across  the  shock  front  are 


pi(ui-c)  =  p2(u2-e) 


(3.4) 


pl(Ui-C)a  +  px  =  p2(u2-C)2  +  p2 


(3.5) 
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S1  =  I  (Pi  +  p2)(£-  •  ^ 


(3.6) 


vhere  subscripts  1,  2  refer  to  conditions  in  front  of  and  behind  the 

shock  front  respectively. 

(b)  Similarity  Assumptions  for  a  Strong  Explosion 

The  approach  used  to  determine  the  self -similar  solution  will 
be  completely  formal,  and  will  employ  a  method  commonly  used  in  the 
solution  of  partial  differential  equations.  That  is,  a  particular 
solution  of  the  partial  differential  equations  is  assumed  reducing  them 
to  ordinary  differential  equations  which  are  then  solved.  This  treat¬ 
ment  will  be  followed  by  a  brief  c cement  on  dimensional  considerations. 
We  seek  a  means  by  which  we  may  obtain  an  analytic  solution  of  the 

problem  of  a  strong  explosion  in  a  continuous  medium.  In  order  to  do 

-a 

so  we  introduce  a  variable  l  =  rt  and  the  following  assumptions  : 


u  =  tP|  U(|) 

p/p  -  tV  P(t)  (3-7) 

P  =  tBD(|) 

where  a,p,e  and  6  are  parameters,  and  U,  D,  P  functions  to  be 
determined.  By  introducing  these  functions  we  shall  show  that  the  par¬ 
tial  differential  equationsof  motion  (Eqs.  (3«l)  -  (3*3))are  reduced  to 
ordinary  differential  equations.  The  variable  |  defines  geometrically 

a  family  of  surfaces  5  =  constant  in  the  r,t-plane,  one  of  which  is 

(18) 

the  shock  front.  (ThiB  has  been  shewn  experimentally  v  ' 


to  be  a  very 


2k 


self -similar  assumptions  are  valid. 

The  notion  of  an  equation  of  state  of  a  medium  implies  thermo¬ 
dynamic  equilibrium  and  can  be  described  by  a  functional  relationship 
between  any  three  of  the  variables  p,  V,  T,  S,  e.  Thus  there 
exists  several  forms  of  the  equation  of  state  of  a  pure  substance. 
Consider  the  Helmholtz  work  function 

F  =  e  -  TS 

From  this  and  the  first  and  second  laws  of  thermodynamics  follows 

dF  =  de  -  TdS  -  SdT  =  -  pdV  -  SdT 

Now  dF  Is  an  exact  differential  of  the  thermodynamic  potential 
function  and  hence 

A 

5v5r 

Since 

(gy^T 

it  follows  that 

de  =  (§|)v  dT  +  (|f)T  dV 

=  Cy  dT  +  [T(^)t  "  PldV 


"  ^5v^t  =  (§r}v 
-  T(|)T  -  P  -  T(fe_  -  P 
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An  ideal  gas  is  defined  as  a  medium  'which  possesses  no  inter- 
molecular  forces  of  restraint;  this  implies  that  the  internal  energy 
is  a  function  of  temperature  alone ;  hence  =  0  and 


p  =  Tg'v)  =  Tf  (p)  (4,1) 


Now  a  solid  possesses  energy  due  to  the  presence  of  intermoleeu- 
lar  forces  and  a  thermal  energy  assoc lated  with  molecular  vibrations. 
Hence  even  at  absolute  zero  the  energy  of  a  solid  is  non-zero.  Thus,, 
we  can  express  the  energy  of  a  solid  mathematically  as 

e  =  f(v)  +  Tg(V)  (4,2) 


Another  formulation  for  the  higher  pressure  range  called  the 

(10) 


Mie-Griineisen 


equation  is 

P  '  P„(p) 


e(p,p)  -  e  (p)  = 


P^"P.) 


(4,5) 


The  subscript  c  refers  to  the  cohesive  contribution  and  f1  is  the 
Griineisen  constant  which  depends  weakly  on  p. 

By  rewriting  equation  (4  5)  in  another  form  the  measured  shock 
wave  data  along  the  Hugoniot  curve  can  be  used  to  calculate  the  co¬ 
hesive  contribution. 


eH(p)  -  ec(p) 


PH(P)  -  PC(P) 

p  r  (p) 


(4  4) 
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could  be  taken  into  account  but  has  been  dropped  for  convenience. 

For  fixed  parameters  a,  P,  €,  6  the  solution  P  =  P(u)  of 
the  governing  differential  equation  (3-20)  is  seen  to  depend  on  y  ■ 

The  value  of  y  determines  the  location  of  the  singular  points  and 
other  significant  aspects  to  be  discussed  in  more  detail  in  the  next 
section. 

In  determining  the  single  value  of  y  to  be  used,  we  must  real¬ 
ize  that  the  y  which  best  describes  the  state  for  one  point  in  the 
medium  will  certainly  be  )  poorer  approximation  for  another  point. 

If  a  straight  line  is  fitted  the  shock  wave  data,  from  Ref.  (19) 
and  is  considered  as  representing  only  the  high  pressure  range,  an 
equation  of  state  for  the  metals  considered  is  obtained.  These  are 
actually  Hugoniot  curves,  but  will  be  used  to  approximate  the  isen- 
tropes  of  the  metals  (See  Huang  <20^). 

We  will  now  consider  a  particular  equation  of  state  p=A( ( p/p^ )  ^  -1 ) 

y 

=  Aft  -  B  called  quasi -polytropic  in  this  thesis  „  This  equation  is 
very  similar  to  the  polytropic  equation  of  state  and  is  used  by  many 
authors  to  describe  the  equation  of  state  for  dense  media,  since  p  ■  0 
does  not  imply  p  =  0  as  for  gases.  Comparison  is  made  between  the 
polytropic  and  quasi -polytrqpic  equations  with  respect  to  relations 
predicted  by  each  equation  or  shock  front  relations,  sound  speed,  and 
thermodynamic  relations.  The  inability  of  the  quasi -polytrqpic  equa¬ 
tion  to  be  used  in  a  self -similar  theory  is  also  exhibited. 

y 

The  sound  speed  for  a  FT  gas  p  *=  Ap  is 

2  y  “X 

c  =  dp/dp  =  yAp' 
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d 


CURVE  FOR  ALUMINUM 
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c  = 


dB 
_ e 

dt 


it  (61 *a> 


or 


C  =  2/5  (1*51.28)  t‘5/5om/i 
=  l8l  t-3/5  cm/i 


sec 


'sec  . 


(4.20) 


The  physical  variables  immediately  behind  the  shock  wave  are 
given  oy  equations  (4.14)  through  (4.l6).  That  is 


"2  ■  ik  c  -  rhi  a8i)f3/5 

=  42.0  t  3//^  cm/sec 


(4.21) 


and 


p2  =  Z75  P1  =  3’518  g/cc 

p2  =  BTF  (2-785)  U8l  t'3/5) 
=  2.13  x  10"5  t'6/5  kb 


(4.22) 


(4.23) 


In  order  to  describe  the  distribution  of  the  variables  through 
the  medium  between  the  shock  wave  and  cavity  surface,  we  turn  to  the 
solution  of  the  differential  equation  (5-20)-  The  solution  of  the 
blast  wave  problem  was  exhibited  as  P  =  P(U).  It  was  shown  that  the 
functions  P  =  P(§),  U  =  U(|),  D=D(|)  could  be  obtained  by  quad¬ 

rature  of  Eqs.  (3*15)>  (3-18)  and  (3*19) •  However,  it  was  found  to  be 
more  convenient  to  proceed  in  a  different  manner. 


VELOCITY  DECAY  ON  CAVITY  SURFACE  FOR  FOUR  METALS 
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radius  to  minimum  values  at  the  cavity  boundary.  We  note  that  the  equa¬ 
tions  of  state  (Fig.  4)  which  have  been  used  for  the  calculations  have 
a  lower  limit  of  about  p  =  100  kb.  This  could  also  be  used  as  a  cri¬ 
terion  for  shock  termination  (point  D,  Fig.  5).  This  point  is  reached 
in  0.7  M  sec.  These  conditions  thus  determine  a  roughly  parallelogram 
shaped  region  ODPE  in  the  r,t -plane  for  the  validity  of  the  assumed 
theory.  Note  that  the  cavity  has  only  expanded  0.1  cm  during  this 
period,  which  is  about  l/30  of  the  total  observed  increase  in  radius. 

We  are  thus  justified  in  referring  to  the  shock  process  as  impulsive, 
i.e.,  the  later  stages  of  the  process  are  insensitive  to  many  features 
of  the  shock  part.  Hence,  the  self -similar  solution  remains  valid  at 
least  for  the  analysis  of  the  shock  zone. 

The  asymptotic  characteristics  of  the  solutions  are  thus  not  of 
direct  physical  interest  since  they  do  not  apply  to  the  problem  much 
beyond  the  region  described  above  The  expansion  zone,  headed  by  a 
wave  travelling  with  the  dilatational  wave  velocity  goes  on  for  at 
least  100  p.  seconds,  during  most  of  which  the  metal  continues  to  move 
by  fluid  or  plastic  flow. 

The  final  cavity  radius  attained  is  of  great  interest  to  the 
general  problem  as  this  value  is  directly  observable  on  the  specimens 
after  blast.  In  principle  the  prediction  of  this  radius  should  afford 
a  test  of  any  theory,  but  the  matter  is  not  so  direct  as  this,  since 
several  theories  are  involved.  It  is  now  evident  that  the  cavity  for¬ 
mation  process  is  complicated.  It  starts  under  one  theory  (in  which 
the  state  of  the  metal  is  fairly  well  established)  but  terminates  in  a 
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ed  elastic  wave  and  the  plastic  loading  or  unloading  wave.  Another 

possible  explanation  is  based  on  assuming  that  the  motion  is  quasi - 

(  2k) 

static  and  fully  plastic.  Hill  '  shows  that  the  displacement  of 
the  inner  surface  increases  with  an  increase  in  outer  radius. 

Figure  15  is  a  sketch  of  a  post-shot  section  of  the  sphere. 

Notice  the  three  regions  marked  on  the  sketch  labeled  A,  B,  C.  Sec¬ 
tion  C  is  a  narrow  region  surrounding  the  cavity  in  which  there  are 
few  ciacks  and  during  deformation  the  metal  appears  to  flow  much  as 
a  fluid.  In  region  B  the  material  exhibited  many  small  shear  and 
tension  failures.  Zone  A  appeared  to  be  dominated  by  reflective  ef¬ 
fects  and  showed  several  "scabbing"  type  failures. 

In  all  of  the  tests,  the  plug  was  blown  out.  It  was  concluded 
that  this  did  not  effect  the  final  cavity  size  significantly  as  the 
sphere  had  to  expand  throughout  to  free  the  plug.  However,  the  pre¬ 
sence  of  the  plug  was  concluded  to  cause  the  asymmetric  deformation 
of  the  sphere. 

The  quartz  disc  crystals  used  to  signal  the  time  of  arrival  of  the 
wave  were  only  able  to  measure  the  small  amplitude  waves.  The  crystals 
would  not  sustain  the  violent  conditions  in  the  shock  zone. 
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CHAPTER  VI 
CLOSURE 

6.1.  Summary  and  Conclusions 

The  behavior  of  shock  waves  in  metals  is  a  subject  of  current 
scientific  and  engineering  interest .  There  has  been  a  considerable 
amount  of  experimental  research  end  a  lesser  amount  of  analytical 
effort  in  this  field.  Seme  of  the  most  successful  approaches  to  the 
problems  of  shock  waves  in  solids  bare  been  based  on  the  so-called 
"hydrodynamic'1  model  0 

In  Chapter  III  of  this  thesis  the  theory  of  self -similar  motion 
is  reviewed  and  discussed  The  equations  for  spherical  motion  of  a 
compressible  in viscid  fluid,  subjected  to  a  strong  shock  are  solved. 

One  of  the  main  contributions  of  this  thesis  is  contained  in 
Chapter  XV  In  this  chapter  the  problem  of  the  dynamic  short  time 
behavior  of  a  ductile  metal  surrounding  a  spherical  detonated  charge 
of  high  explosive s  is  solved  for  the  first  time.  The  problem  is  ap¬ 
proached  using  the  theory  of  self -similar  motion  presented  in  Chapter 
III  The  trajectories  of  the  shock  waves  and  cavity  surface  are  ob¬ 
tained  as  well  as  values  for  the  variation  with  time  of  pressure  and 
velocity  on  the  cavity  surface.  Cavity  formation  in  other  metals  is 
discussed  and  values  for  shock  and  cavity  surface  trajectories  are  ob¬ 
tained.  Pressure  and  velocity  variation  with  time  on  the  cavity  sur¬ 
face  are  also  presented. 

A  critical  analysis  of  the  theory  of  self -similar  motion  as  it 
applies  to  shock  waves  in  solids  is  presented.  The  restrictions  im- 
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posed  on  its  use  by  equations  of  state  and  assumptions  of  constant 
shock,  strength  are  discussed.  An  equation  of  state  called  quasi-poly- 
tropic  is  considered  and  its  incompatibility  with  the  theory  of  self¬ 
similar  motion  exhibited. 

A  method  which  accounts  for  one  of  the  non-similar  aspects  en¬ 
countered  at  large  times  for  spherical  shock  waves  in  solids  is  dis¬ 
cussed  and  its  application  to  the  cavity  expansion  problem  outlined. 

The  results  of  the  analysis  for  aluminum  show  that  the  "shock 
zone"  or  range  of  validity  of  the  self -similar  solution  is  very  small. 
The  reasons  for  this  are  (l)  after  only  0-5  microseconds  the  shock 
wave  velocity  predicted  by  the  self -similar  solution  has  fallen  to  that 
of  the  elastic  wave  velocity  in  the  material,  (2)  in  only  0.7  micro¬ 
seconds  the  pressure  behind  the  shock  front  has  fallen  to  100  kb  (which 
has  been  chosen  as  the  lower  bound  for  the  range  of  validity  of  the 
equation  of  state  of  aluminum).  At  this  time  the  cavity  has  expanded 
0.1  cm  which  is  only  l/30  of  the  total  increase  in  radius. 

Thus  the  material  behavior  is  referred  to  as  impulsive,  that  is, 
the  latter  stages  of  the  process  are  insensitive  to  many  of  the  fea¬ 
tures  of  the  shock  part. 

The  final  cavity  radius  attained  is  of  great  interest  to  the  gen¬ 
eral  problem  as  this  value  is  directly  observable  on  the  specimen  after 
blast.  In  principle,  the  prediction  of  this  radius  should  afford  a  test 
of  any  theory,  but  the  matter  is  not  so  direct  as  this,  since  several 
theories  are  involved.  It  is  now  evident  that  the  cavity  formation  pro¬ 
cess  is  complicated.  It  starts  under  one  theory,  in  which  the  state  of 
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static  and  the  sphere  completely  plastic.  Under  these  assumptions  Hill 
shows  that  the  inner  radius  increases  with  the  outer  radius.  Other 
things  remaining  equal. 

The  damage  suffered  hy  the  spheres  is  seen  to  he  rather  sharply 
delineated  into  3  categories.  The  first  is  the  zone  immediately  next 
to  the  cavity.  In  this  zone  there  were  few  shear  cracks  and  the  ma¬ 
terial  appeared  to  have  been  deformed  in  a  fluid-like  manner. 

In  the  second  zone  this  was  not  the  case.  There  were  numerous 
shear  and  tension  cracks  and  damage  was  severe. 

The  outer  zone  was  dominated  by  reflective  effects.  Most  of  the 
failures  in  this  zone  were  "scabbing"  type  failures 

.Another  interesting  fact  determined  from  the  experiments  is  that 
there  is  definitely  seme  oscillation  of  the  sphere. 

In  all  the  tests  the  plug  was  blown  out.  It  was  decided  that  this 
did  not  affect  the  cavity  size  to  a  great  extent,  since  it  only  occurr¬ 
ed  after  the  sphere  had  expanded  sufficiently  to  allow  its  release. 

6.2  Suggestions  for  Further  Research 

It  Is  obvious  In  considering  the  analysis  of  later  stages  of  the 
cavity  expansion  process  that  a  criterion  is  needed  for  marking  the  end 
of  the  self -similar  solution- 

Several  such  mechanisms  have  been  suggested  such  as  temperature, 

(10)  (Q) 

material  strength  '  and  energy  but  none  are  completely  satis¬ 
factory. 

To  obtain  a  better  approximation  of  the  actual  material  behavior 
for  later  time,  one  must  account  for  the  non-similar  aspects  of  the 
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APPENDIX  A 

SINGULAR  POINTS  OF  ORDINARY  DIFFERENTIAL  EQUATIONS 

The  singular  points  of  an  ordinary  differential  equation  dv/dx  = 
f(x,v)/g(x,v)  are  defined  to  "be  those  points  for  which  dv/dx  =  ^  . 

It  is  necessary  to  study  carefully  the  singular  points  of  the  dif¬ 
ferential  equation  (Eq.  (3  20)  ),  ~  =  >  111  order  to 

construct  solution  curves  using  numerical  methods. 

Poincs.re  shewed  that  the  differential  equation 

dv  _  ax  +  bv'  *  (A.l) 

dx  cx  +  dv  +  g2(x,d) 

in  which  the  constants  a,  b,  c,  d  are  such  that  A  =  ad  -  be  /  0 

2  2 

and  in  which  fg  and  g2  vanish  like  x  +  v  as  x,v  -»  0,  has  as 
its  only  singularities  those  of  the  much  simpler  equation 


dv  '  ax  -t-  bv 
dx  ~  cx  +  dv 


(A.2) 


It  is  obvious  that  these  occur  only  at  the  origin  rx,v)  =  (0,0). 
However,  in  order  to  study  a  singular  point  (^0>  PQ)  a  change  of  var¬ 
iables  can  be  made  by  letting  x  =  U  -  UQ,  v  =  P  -  Pq.  Thus  in  the 
neighborhood  of  an  isolated  singularity (UQ,?o)  the  differential  equation 
(3-20)  can  be  written  .  The  behavior  of  the  solution 

curves  in  the  neighborhood  of  the  point  is  then  determined  by  the  val¬ 
ues  of  a,  b,  c,  d. 


77 


For  aluminum  (7  «=  7-6)  the  point  (S4)  is  given  by  U  =  0.0916, 
P  =  0.02884.  The  remaining  roots  are  imaginary. 

A  summary  of  the  results  of  the  classification  of  the  nodes  is  as 
follows : 


.nguJ 

Lar  Point 

Classification 

SI 

(0,0) 

stable  node 

S2 

(a,o) 

saddle  point 

S3 

(1,0) 

stable  node 

S4 

(.0915,  .02884) 

unstable  node 

From  numerical  calculations  it  appeared  that  the  solution  curve 
for  7  *>7-6  should  end  at  S2.  However,  this  was  a  saddle  point  of 
the  differential  equation  and  curves  starting  near  S4  never  went  into 
S2.  (Theoretically  we  know  there  is  only  one  curve  which  does  so  and 
this  cannot  be  found  numerically.)  In  order  to  obtain  a  better  approxi¬ 
mation  to  a  solution  curve  joining  the  points  S2  and  S4  the  following 
procedure  was  used. 

Near  U  =  a  we  have 

PcPo+  <§>U=a  <U  -  a)  +  ••• 

Dividing  both  numerator  and  denominator  in  the  differential  equation 
(A. 3)  by  U  -  a  and  putting 

m  =  P/U-a  (A. 5) 


we  may  write 
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APPENDIX  B 
NUMERICAL  PROGRAM 

The  problem  programmed  for  the  computer  was  to  obtain  the  traject¬ 
ory  of  the  shock  front  and  particle  paths  as  well  as  distribution  of 
physical  -variables  from  the  differential  equation  and  a  starting  point, 
(see  Appendix  A)  together  with  boundary  conditions  and  kinematic  rela¬ 
tions  . 

In  order  to  avoid  the  difficulties  caused  by  vertical  or  horizon¬ 
tal  tangents,  the  differential  equation  (A. 5)  was  written  in  the  follow¬ 
ing  parametric  form, 

—  .  —  (R  +  PS)  -  dU  (B.l) 

y  P2(NfPQ)B  +  (R+FS)2 

. . —  P(N+PQ)  =  dP  (B-2) 

y  P2(N+PQ)2  +  (R+PS)2 

where 

ds  =  y  dP2  +  dU2  =  y  1  +  (||)2  dU  (B.3) 

Expressions  (B.l)  and  (B.2)  are  well  behaved  except  at  S4  and  S2, 
which  are  carefully  avoided. 

The  accuracy  of  the  calculation  is  seen  to  depend  on  the  choice  of 
ds  (other  things  remaining  fixed).  However,  ds  has  no  physical  mean- 
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TABLE  4 

Program  Symbol  Table 


Program  Symbol  Our  Notation  or  Name 


Mode 


A 

a 

FLT 

B 

P 

tt 

C 

6 

It 

G 

7 

If 

U 

U 

If 

P 

P 

If 

DS 

ds 

II 

X 

6 

II 

D 

P 

II 

IMAX 

num.  of  iterations 

FIX 

EP 

potential  energy 

FLT 

EK 

kinetic  energy 

II 

E 

E(t) 

II 

FN 

N(U) 

II 

FQ 

Q(U) 

If 

FR 

R(U) 

II 

FS 

s(u) 

II 

FT 

P(N+PQ) 

II 

FW 

R+PS 

II 

FX 

ds/y  P2(NfPQ)2+(R+PS)2 

II 

DP 

dP 

It 

DU 

dU 

It 

FY 

(u-a)7-rP 

If 

DX 

If 

DD 

dp 

II 

PD 

pp 

II 

DEP 

II 

DEK 

II 

DO 

Pi 

II 

R 

R 

II 

St.  47 

initial  U 

tl 

St. 48 

initial  P 

II 

St.  54 

initial  p 

II 

EX 

1/a 

II 

T 

t 

II 

V 

u 

II 

Pres 

P 

It 

DT 

dt 

II 
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spherical  charge  of  high  explosive,  subsequent  to  the  detonation  of  the 
charge . 

The  problem  is  approached  using  a  similarity  theory  (a  one  dimen¬ 
sional  problem  in  vhich  the  motion  can  be  described  by  one  independent 
variable  and  hence  the  partial  differential  equations  of  motion  are  re¬ 
placed  by  ordinary  differential  equations . )  A  more  complete  explana¬ 
tion  of  this  type  of  theory  is  given  in  reference  (4)  page  146.  Values 
for  the  velocity,  pressure  and  density  of  the  material  are  found  for 
the  early  stages  of  the  process  The  results  are  compared  with  those 
of  experiments.  In  addition,  a  critical  study  is  made  of  the  range  of 
validity  of  similarity  type  solutions  as  they  apply  to  shock  waves  in 
solids .  A  method  of  extending  the  self -similar  theory  is  described 

Equations  of  state  of  solids  are  discussed  and  their  compatibil¬ 
ity  with  the  theory  of  self -similar  motion  studied.  A  particular  type 
equation  of  state  called  quasi -polytropic  in  this  work  is  analyzed  as 
to  its  compatibility  with  the  similarity  solution.  Comparison  is  made 
between  the  shock  relations  obtained  using  this  quasi -polytropic  rela¬ 
tion  and  the  polytropic  equation  of  stare 

Experiments  on  thick  aluminum  spheres  are  described  and  the  re¬ 


sults  are  discussed. 
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Several  interesting  approximate  solutions  for  large  elastic -plastic 
deformation  are  also  discussed  as  well  as  a  numerical  approach  using 
the  theory  of  characteristics. 

2.2.  Similarity  Solutions  of  Equations  of  Gas  Dynamics 

Many  articles  have  appeared  in  the  literature  on  similarity  solu¬ 
tions  of  the  equations  of  gas  dynamics.  The  particular  equations  which 
appear  in  this  thesis  were  first  presented  by  G.  I.  Taylor 

Taylor  worked  out  a  few  numerical  solutions  as  well  as  seme  ap¬ 
proximate  analytical  solutions  for  an  explosion  in  a  perfect  gas  with 
7  ranging  from  1.2  to  I.67.  Subsequently  or  simultaneously  analytic 

solutions  appeared  by  Sedov  J.  L.  Taylor  Latter  and 

(7) 

Sakural  ' 1 ' .  Sedov  presents  the  most  complete  discussion  of  the  solu¬ 
tion  curve  in  the  plane  of  the  non-dimensional  pressure  and  velocity 
functions.  He  discusses  the  singular  points  and  their  significance  in 
the  physical  problem. 

/  a) 

However,  it  is  amazing  that  Courant  and  Friedrichs'  ,  who  ob¬ 
viously  were  not  aware  that  the  exact  solution  existed,  were  able  to 
deduce  nearly  as  complete  a  description  of  the  field  of  solution  curves 
as  Sedov. 

2.3.  Self -Similar  Solutions  of  the  Problem  of  Cratering  Due  to  Hyper¬ 
velocity  Impact 

As  was  mentioned  earlier,  the  problem  of  blast  waves  in  spheres 


and  cratering  due  to  hypervelocity  impact  have  much  in  common.  Conse¬ 
quently,  the  approach  used  in  this  thesis  is  similar  to  that  used  by 
several  investigators  of  hypervelocity  impact  phenomena. 
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good  assumption. 

We  new  explore  these  solutions  mathematically  hy  substituting  the 
expressions  (3. 7)  into  the  equations  of  motion,  giving  respectively. 


|tP_1  [pu  -  a  (|U'+U)  +  tP_a+1  U(|U'+U)  +  t6-01"^1 
(2P  +  i  PD'/D  +  |P')]  =  0 


(3.8) 


C  n 

t  [(U-a)  £D’  +  BD  +  ( |U'  +  3U)D]  =  0 


(3.9) 


t&4-€-l|-_  ^2  p(BD^D.)  +  ^2  jjp 

-  a(2DP+|D'P  +  £DP')1  +  tP-a+1[  -  7|2PD'  (3.10) 

+  |(2DP  +  |D'P  +  £DP ' ) ]  |U]  =  0 

7 

where  the  equation  of  state  p  =  Ap  has  been  used.  This  will  be  dis¬ 
cussed  in  detail  later  in  the  thesis.  These  equations  are  new  in  a 
form  where  it  is  possible  to  eliminate  the  explicit  factor  t  by  prop¬ 
erly  choosing  the  exponents,  thereby  leaving  a  system  of  functions  of 
one  independent  variable  |.  This  is  accomplished  by  letting 


s  o  that  t 


e  =  2p 
P-0M-1 


=  t 


p  =  a-l 
e-a-^fl 


Dividing  by  |t 


P-1 


*=  1* 
J3-1 


5+e-l 


(3-11) 


(which  are  not  zero  for  t>  0) 


and  simplifying  the  resulting  equations,  we  obtain. 


E(t)  =  4*  f  1  (t**6  FI  +  |  t&+2P  tJ2)62D(|)r2dr 
r 

o 

Substituting  r=|ta ,  dr«tad|  into  the  preceding  equation  and  noting 
that  t  is  constant, 

E(t)  =  bn  t8+5a’2  f  (  JL  +  •  (3.24) 

gQ  7-1 

Since  the  integral  is  independent  of  t,  we  make  the  energy  independ¬ 
ent  of  time  by  satisfying  the  relation 

6+50-2  =0  or  6=2-50  (3-25) 

We  shall  narrow  down  the  number  of  parameters  by  examining  the 
compatibility  of  our  solution  with  the  basic  Rankine-Hugoniot  con¬ 
ditions  across  a  shock  front.  If  the  undisturbed  and  disturbed  medium 
parameters  are  u^,  p^,  p^  and  Ug,  Pg,  p2,  respectively  and  the 

shock  wave  velocity  is  C,  then  the  relations  for  (see  (8)  p.123-4), 
conservation  of  mass,  momentum  and  energy  across  the  shock  front  are. 


P1(u1-C)  =  pg(Ug-C)  =  m 


(3.26) 


PgUg^g-C)  -  P1U1(U-L"C)  =  P-^Pg  (3.27) 

P2(|u22+c  )  (u2-C)  -pi(|ui2+c)  (ui-C )  =  P^-pgUg  (3.28) 

where  e  =  p/p  for  a  polytrqpic  medium. 

When  the  undisturbed  state  is  a  medium  at  rest,  with  u^O,  these 


equations  reduce  to 


20 


plete  set  of  exponents  is  new 

a  =  2/5  €  -  -  6/5 

P  -  -  3/5  6=0 


(3  42) 


If  ve  new  turn  to  purely  dimensional  considerations  the  following 
are  recognized  as  the  system  of  fundamental  parameters  influencing  the 
motion  of  the  medium:  p^,  p^,  Eq,  r,  t,  y.  From  these  we  can 

form  only  three  independent  dimensionless  quantities :  7  ; 


1  r 


3/6 1 


E 


r/5'r~275~'  *  T  =  "1/3 - \J2~  ”  For  a  8tron8  shock  (p1«p2) 

1  Eo  P1 


the  system  reduces  to  only  7  }  X  "E,^3  /p~^3  =  rt  2^.  If  we 

take  7  to  he  constant,  we  see  that  the  motion  can  he  described  by 
only  one  independent  variable,  X  ,  which  agrees  with  our  earlier  find¬ 
ings  . 

We  can  then  write  (see  Sedov,  Ref.  (4)  ) 


u  =  (r/t)  V(l)  ,  p  =  tS_H2  jt(x) 

r 


=  tS  R(X)  a  =  [MLV] 


(3=43) 


This  system  of  equations  differs  frcm  the  ones  previously  derived 
only  in  the  dimensional  sense.  Here  it,  V  and  R  are  all  dimension¬ 
less  functions  of  X  and  are  related  to  the  functions  P,  U,  D  in 
the  following  way:  U  =  V,  R  =  D/p.^  P  =  "P1AI  “  |  • 
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Subtracting  (4.4)  from  (4.5)  we  obtain 


e  -  e  (p)  «  p  PH^ 

H  Trl  P) 


The  Mie-Griineisen  equation  can  be  rewritten  in  the  form 


e 


P 

pT’TpT 


-  A(p) 


and  we  find 


A(p) 


Pc(p)  Pjj(p) 

p  r  (p)  '  ec(p)  =  pp(f>)  ‘  eH  p 


(4.5) 


The  last  expression  can  be  found  from  shock  wave  data  and  in  this  way 
the  magnitude  of  A(p)  can  be  determined. 

It  has  been  found  that  the  only  form  of  the  equation  of  state 
compatible  with  the  self -similar  flow  assumption  i3  of  the  form 


e  =  P<p(p)  • 

We  see  that  the  leading  term  of  the  Mie-Griineisen  equation  is 
p/p P(p)  =  pcp(p)  and  we  are  forced  to  use  only  this  part  of  the  equa¬ 
tion  in  the  self -similar  solution.  Thus  the  assumption  is  accurate  so 
long  as  the  term  A(p)  is  small  compared  to  p//-7(p)-  It  should  be 
pointed  out  that  it  does  not  follow  that  the  fluid  assumption  is  in¬ 
valid  when  A(p)  is  not  small  compared  to  v/pP  (p)* 

The  Griineisen  factor  P  has  been  assumed  to  be  a  constant  equal 
to  (7-1).  This  is  not  absolutely  necessary  since  the  variation  of  P 
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when  the  values  (v^,p^)  on  the  other  side  are  prescribed.  For  a  poly- 
tropic  equation  of  state  it  is  a  rectangular  hyperbola  (Fig.  l). 

For  a  QPT  medium 

H(V,p)  -  (p+7B)V  -  ~p—  (p.  +  7B)V 

2pT  2p 

pi  +  P\ 

+  (v-vx)  ' > 


2^  =  (l-p2)(pV)  +  (l-p2)7BV  -  (l-p2)?^ 


o  o  P  +  P-, 

-  (1-p  )7BV1  +  2u^(V-V1)(— g-  ~) 


2v2n  =  (v-n^Jp  -  (v1-p2v) 


+  (1-p  )  rBiv-v^ 


(4.8) 


Figure  (l)  shows  the  Hugoniot  for  FT  and  a  QPT  medium  plotted  for  the 
purpose  of  comparison  (p1  has  been  assumed  equal  to  zero). 

For  the  case  of  a  shock  moving  into  a  gas  at  rest  we  obtain 

(V-p^tp  +  (l-p2)7B]  =  v1  (i-p2)b. 


This  can  be  rewritten  as 


tt  ?7.B  ( v  -v  ) 

2  ' V2  1; 

-  P  V-, 


2  ^  1 
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FIG.  I  -  HUG0NI0TS  FOR  A  POLYTROPIC  AND 
QUASI-POLY  TROPIC  MEDIUM 
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Following  the  same  procedure  for  Pfl^)  and  we  find 


P(|1)  =  2  a2  p2  /  (7+1) 
D(ll)  =  Px  /  42 


(4,18) 


For  aluminum  with  y  -  7  6  this  supplies  an  initial  point  in  the 
P-U  plane.  The  solution  curve  is  the  unique  integral  of  (5.20)  which 
passes  through  this  point  and  the  singular  point  D,  as  shown  in  Fig. 
3. 

We  need  next  to  prescribe  the  initial  conditions  at  the  boundary 
of  the  cavity  and  the  total  energy  supplied  by  the  explosive.  We  will 
see  how  these  allow  us  to  determine  1^,  and  consequently,  the  motion 
of  the  shock  front. 

From  private  communication  with  R.  E  Shear  of  Ballistic  Research 
Laboratories,  Aberdeen  Proving  Grounds,  we  found  that  for  pentolite 
under  confinement  the  detonation  pressure  would  he  around  400  kb  For 
convenience,  a  value  of  pressure  =  384  kb  was  used  in  this  work.  This 
pressure  predicts  the  proper  density  ratio  using  the  Hugonlot  equation 
(Fig  4)  as  the  equation  of  state  for  the  material  at  the  cavity  sur¬ 
face. 

The  next  step  is  to  assure  that  when  the  shock  wave  is  at  the  cav¬ 
ity  surface,  the  pressure  and  density  have  the  right  values,  namely 

2 

P  =  (^)  P(lx)  »U±) 
o 


(4.19) 
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In  this  expression  p  p,  r  ,  P  and  D  are  known,  and  t  is  de- 

>  o  o 

termined,  i.e., 

-  r0  (p/p  P^l>  1/2 

The  values  used  to  calculate  t  were 

o 

Px  =  2.785  g/cc 
p  =  384  kb 
rQ  =  I.698  cm 
P^)  «  0.0286 

D^)  =  3.5182  g/cc 

as  indicated  in  Fig.  11  and  Table  2.  Using  these  values  we  find  tQ  = 
0.868  p  sec. 

This  value  of  t  is  the  "delay"  time  to  assure  the  proper  position 
of  the  shock  front  at  the  initiation  of  the  expansion  process. 

With  the  values  of  r  ,  t  known,  is  calculated  from  the 

definitive  relation  |  «=  rt~a.  Then 

=  r  t  =  451.28  cm/see2/^ 

’loo  ' 

I 

The  velocity  of  the  shock  wave  is  given  explicitly  by  the  ex¬ 


pression 


4l 


Frofn  Table  1  ve  see  that  the  energy  does  tend  to  be  a  constant  and  we 
have  Just  shown  that  this  limit  is  the  value  of  the  energy  integral 
(3.25). 

If  we  now  suppose  that  all  (or  any  known  fractional  part)  of  the 
energy  given  up  by  the  explosive  is  transmitted  into  the  solid,  then 
the  shock  process  could  be  terminated  when  the  energy  reaches  the 
amount  available  It  Is  not  possible  to  determine  a  precise  point  of 
time  because  of  the  asymptotic  way  in  which  the  energy  increases.  It 
is  seen,  however,  that  E  reaches  90 $  of  its  ultimate  value  in  5  micro 
seconds,  which  is  a  very  short  time  compared  with  the  expansion  process 

The  post -shock  expansion  presumably  must  take  place  under  con¬ 
stant  energy  conditions  for  a  "long"  period  of  time,  until  it  is  dissi¬ 


pated  by  viscosity  of  the  flow,  elastic  waves  and  other  side  effects. 


For  energy  available  in  the  explosive.  Shear  '  J  gives  the  value 


■Wi  ■ 1-152  k  “1/e 

which,  for  our  explosive  weight  of  0.07  It#  gives 
*rotoi  ■  126  *  lo12  “«*• 


Our  calculated  asymptotic  value  (Table  l)  from  the  energy  integral 
coraes  to  96$  of  this  value 
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u/u2  =  l/6x  -  r/k0  =  p/p2 

p/p2  -  (l/^)3 


(4.24) 


A  value  of  7=7  serves  to  separate  the  domains  of  the  function 
U.  For  7  <  7  i 

w  * 0  s  irfer 

and  for 

7  >  7  <  5T7HT  f  D  -  2/5  ' 

The  value  U  -  corresponds  to  the  shock  wave  in  a  gas 

due  to  a  point  explosion  and  2/5  7  corresponds  to  the  center  of  the 
explosion.  However,  in  the  cavity  expansion  it  is  meaningless  and 
hence,  lower  hound  of  the  domain  of  U  would  he  some  number  greater 
than  2/5  7* 

For  7  >  7  the  domain  of  the  function  U  is  bounded  above  by 
U  =  2/5  which  corresponds  to  the  edge  of  a  so-called  vacuum  for  a 
point  explosion  in  a  gas  and  below  by  the  shock  front.  The  value  of 
U  =  2/5  would  not  necessarily  apply  to  the  cavity  problem.  Fig.  1 6 
shows  the  behavior  of  the  solution  curves  for  variable  7* 

However,  it  is  significant  that  the  solution  curve  P  =  P(U), 
which  describes  the  process,  has  entirely  different  behavior  in  the 
two  ranges.  (See  Figs.  10,  16).  Figure  10  shows  the  changeover  from 
solution  curves  which  go  to  point  B  from  those  which  go  to  D,.  (Fig. 3). 


It  is  fitting  at  this  time  to  mention  that  the  careful  considera¬ 
tion  of  the  singular  points  and  their  physical  meaning  came  about  frcan 
having  the  analytical  solution  of  the  strong  blast  wave  problem.  It 
is  doubtful  if  this  could  have  been  accomplished  by  numerical  computa¬ 
tion  alone.  This  also  applies  to  establishing  the  exact  domains  of 
the  function  U.  However,  it  is  possible  to  determine  the  type  of 
singular  point  that  exists  as  well  as  a  good  description  of  the  gen¬ 
eral  behavior  of  the  solution  curves  without  the  analytical  solution 
of  the  differential  equation. 

4-  5  -  Consideration  of  the  Non -Similar  Aspects  of  the  Blast  Wave 
Problem 

It  is  well  known  that  even  for  the  case  of  spherical  blast  waves 

in  air  that  for  large  time  the  flow  becomes  non-similar.  That  is,  the 

motion  can  no  longer  be  described  by  the  self -similar  assumptions. 

This  is  particularly  true  for  the  case  of  shock  waves  in  metals 

where  use  has  been  made  of  the  strong  shock  approximation.  For  large 

values  of  time  the  density  ratio  Pg/Pp  is  no^  constant  and  the  shock 

2/5 

front  is  no  longer  described  by  R  =  |7t  '  .  In  order  to  extend  the 

S  X 

domain  t  of  the  self-similar  solution  it  is  necessary  to  account  for 
these  non-similar  effects. 

An  approximate  theory  developed  by  Rae  for  the  cratering 

problem  will  be  considered  and  its  use  in  the  cavity  expansion  problem 
described.  The  assumption  is  that  for  any  shock  speed  C  *  Rs  the 
distribution  of  velocity,  pressure,  density,  etc.,  are  the  same  as  for 
the  self -similar  polytropic  equation  of  state  distributions  which  would 


RELATIVE  POSITION  OF  S4  AND  SHOCK  FRONT  FOR 
VARIABLE  7 
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Equation  (4. 30)  describes  the  shock  radius  in  terms  of  shock  speed 
and  is  a  very  useful  relationship.  By  considering  E  as  a  free  para¬ 
meter  Eq.  (4.30)  could  be  used  to  fit  test  results.  On  the  other  hand, 
the  expression  could  be  solved  for  E  and  estimates  could  be  made  of 
the  percent  of  available  energy  present  in  the  solid. 

If  equation  (4.30)  is  solved  for  R  =  R  (R  )  and  the  identity 

s  s  s 

R 

pS  dR 

t  =  J  +  k  (4.31) 

R  8 
s 

is  used,  we  obtain  an  expression  for  Rg  =  Rg(t).  The  value  of  k  is 
determined  fron  consideration  of  the  "delay  time"  discussed  in  Section 

4.3. 

With  the  relations  R  =  R  (t)  and  R  =  R  (r  )  at  our  disposal 

s  s  s  s  s 

it  is  possible  to  obtain  the  distributions  of  velocity  pressure  and  den¬ 
sity  throughout  the  solid.  The  cavity  radius  versus  time  curve  can 
also  be  found.  However,  the  procedure  used  in  Section  4.3  would  have 
to  be  modified. 

We  recall  that  in  solving  the  problem  for  constant  7  that  we  con¬ 
structed  one  solution  curve  P  =  P(U)  which  described  the  functions 
for  all  g.  In  the  ca3e  when  7  is  allowed  to  vary  continuously,  an 
infinite  number  of  solution  curves  must  be  constructed. 
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The  obvious  way  to  avoid  this  is  by  considering  small  increments 
of  time  At  at  the  end  of  which  the  shock  speed  is  calculated.  This 
value  of  Rg  is  used  to  calculate  the  corresponding  7  frctn  Eq. 

(4  25)  and  for  this  instant  the  distribution  of  velocity,  pressure  and 
density  as  well  as  the  particle  velocity  on  the  cavity  surface.  This 
would  be  continued  to  the  limit  of  validity  of  the  equation  of  state 
or  by  using  some  other  criterion  which  terminated  the  solution. 

4.6.  Summary  and  Conclusions 

In  this  chapter  we  have  attempted  to  study  the  cavity  formation 
process  in  several  metals  by  determining  how  the  important  physical 
variables  of  cavity  radius,  velocity,  pressure,  and  density  vary  with 
time  and  position  near  the  cavity.  The  roost  prominent  general  feature 
of  the  whole  process  is  the  short  time  of  the  "shock"  regime  as  compar¬ 
ed  with  the  total  time  of  the  expansion.  One  general  criterion  for 
the  end  of  the  shock  is  furnished  from  the  time  when  the  super-sonic 
velocity  of  the  shock  front  drops  to  sonic  i  e.,  at  the  point  P  in 
Fig.  5 ,  where  the  slope  attains  the  value  for  elastic  disturbances  in 
the  material.  There  is  some  concern  whether  the  self -similar  solution 
can  be  used  beyond  this  point  since  it  would  predict  a  subsonic  shock 
velocity.  However,  it  could  mean  that  a  large  plastic  wave  is  devel¬ 
oped  which  travels  at  a  speed  less  than  the  elastic  wave  velocity. 

This  transition  region  is  highly  uncertain.  This  situation  has  occurr¬ 
ed  after  0.5  microsec. 

We  note  that  the  highest  pressures  and  densities  in  the  metal  are 
located  Just  behind  the  shock  front,  and  trail  off  with  decreasing 
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CHAPTER  V 

EXPERIMENTAL  PROCEDURE  AND  RESULTS 
5-1-  Description  of  Specimens  and  Experimental  Arrangement 

The  experiments  reported  on  in  this  chapter  were  performed  at  the 
Ballistic  Research  Laboratories  of  Aberdeen  Proving  Grounds,  Aberdeen, 
Maryland  The  writer  acted  as  an  advisor  on  the  test  program  and  par¬ 
ticipated  to  a  minor  extent  in  the  testing.  Same  of  these  results  will 

(23) 

appear  in  a  joint  paper  with  a  member  of  the  Ballistic  Research 

Laboratories  staff  to  be  published  later. 

The  specimens  chosen  for  the  tests  were  spheres  made  of  2k  ST  and 
6061 -T4  aluminum  alloys.  The  first  series  of  tests  were  made  using 
spheres  whose  outer  diameter  was  nominally  eighteen  centimeters.  The 
second  series  was  made  using  spheres  of  25  cm  outer  diameter.  The 
inner  cavity  in  both  cases  was  J.k  centimeters.  (See  Fig.  ll).  The 
exact  dimensions  are  given  in  Table  2.  A  threaded  well  was  machined 
in  the  spheres  in  which  a  threaded  plug  with  a  hemispherical  cup  ma¬ 
chined  in  the  end  was  inserted.  This  arrangement  offered  essentially  a 
condition  of  an  isolated  cavity  inside  a  metal  sphere.  A  small  diameter 
hole  drilled  through  the  axis  of  the  plug  provided  a  means  for  running 
the  firing  line  and  ionization  probe  to  the  cavity.  The  end  opposite 
the  hemispherical  cup  protruded  from  the  sphere  and  served  as  a  means 
of  mounting  the  specimen. 

The  explosive  used  was  a  32  gram  sphere  of  Pentolite  which  fitted 
snugly  into  the  cavity  in  the  aluminum  sphere . 


RESPONSE  OF  CAPACITANCE  GAGE  FOR  SMALLER 
SPHERE 


Zone  A  —  No  damage,  other  than  that  of 
reflected  wave 


Zone  B  -  Heavily  damaged 

Zone  C  —  "Fluid"  zone 

Note  —  Dimensions  are  approximate 


FIG.  I  5 


SMALLER  SPHERE  -  POST  SHOT. 


BEHAVIOR  OF  SOLUTION  CURVE  FOR  VARIABLE 
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The  singular  points  of  Eq.  (3.20)  will  now  be  determined  and  class- 

(25) 

ified  according  to  the  notation  and  criterion  of  Stoker  '  .  Revrit- 

ting  Eq.  (3*20) 


dP  P[N(U)  +  PQ(U 
dU  “  R(U)  +  PS(U) 


By  inspection  of  Eq.  (A. 3)  we  find 


(A.3) 


(51)  U  =  0  ,  P  =  0 

(52)  U  =  a  ,  P  =  0 

(53)  U  =  1  ,  P  -  0 


The  location  of  (Sk)  is  determined  by  the  condition 


This  statement  leads  to  the  cubic  equation 
Aj  U5  +  AgU2  +  +  Aq  =  0 

where 

A?  =  27  (1  -  37) 

A2  =  372  (37  -  1)  +  7(7  -  5a)  -  7(38  +  hp)  +  5  +  20 
A1  =  7  (3Q8  +  6otp  -20  -  l*a)  +  (2  -  a)  (8  +  2p) 

Aq  =  -  2a  (8  +  20) 
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TABLE  5  -  FORTRAN  PROGRAM 

COMPILE  RUN  FORTRAN 
C  R-T  CALCULATION  PROGRAM 

100  DIMENSION  I  DENT  ( 2 ) 

1  READ  8OO.A.B.C.G.D0* 1DENT 

800  FORMAT ( 5F10.0.2A5 > 

2  READ  801.X.R.DD.IMAX.PRES 

801  FORMAT < 3F10.0.  I  1Q.1F10.0I 

3  IFUMAX)  44,44»4 

4  PRINT  802 

802  FORMAT  (  1H1  .40X • 32HBRL-PA. STATE  PROJECT  BLAST  WAVES/ 

11H0.48X.15HR-T  CALCULATION/ 

HHOf  15X  •  10HINPUT  DATA) 

5  PRINT  803  .  IDENT .G.R.DO.PRES. IMAX.DD 

803  FORMAT ( 1H0.10X . 1 1HMATER  IAL  «  .2A5/ 

11H0.10X .8HGAMMA  *  *1F5.2/ 

11H0. 10X.17HINIT IAL  RADIUS  *  .1F5.3.3H  CM/ 

11H0.10X.18HINIT IAL  DENSITY  -  .1F5.3.6H  GM/CC/ 

11H0.10X.19HINIT IAL  PRESSURE  -  .1F5.1.9H  KILOBARS/ 

1 1H0 .10X .20HN0.  OF  ITERATIONS  *  .114/ 

11H0. 10X .20HDENS I TY  INCREMENT  *  .1F6.3.6H  GM/CC// 

11H0.1X.1HI  ♦  7  X  ,  1  H  T  .11  X.1MR  .  12X  •  1HV.12X.7HPRES  KB  .  9X  .  1  HD  .  1 IX  .  1HU  •  1  IX 
1.1HP.9X. 1HX .11X.1HE//) 

6  U»A*( 1 • 0— ( G-l #0 )  / tG+1 • 0  )  ) 

7  P*  (  a*A*  (G— 1.0)  /  (  G+l  *0  )  )*(l.O-IG-1.0)/(GU.O)) 

8  D=(G+1.0)/(G-1.0)*00 

9  EP  =  0 

10  EK«0 

11  E  *0 

12  DO  42  1=1  .IMAX 

13  FN»G*U* ( 3 • 0*A-1 .0-2.0*U  )  +  <  3 . 0-A ) *U-2 . 0*A 

14  FQ*<  2.0*B-<G-1.0  )*C)/<U-A)-*-2.0*G 

15  FR*U*  (U  —  A  )  *  (  1.0*"U  ) 

16  FS«C+2.0*B+3.0*U*G 

17  ft=P* ( FN+ P*FQ) 

18  FW-FR+P*FS 

19  F  Y  =  (  U— A  )  *  (  U-A  )  —  G*P 

20  F  X  *-DD*  (  U  —  A  )  /  (  D  *  (  FW+3.0*U*FY)  ) 

21  DP=FT*FX 

22  DU*FW*FX 

23  EX-l.O/A 

24  T“(R/X)**EX 

25  DX  =  X  *  F  Y  *F  X 

26  VCR*  U/T 

2  7  PRES  =  (  R  /  T  )  **2*P*D*(  10 . 0** ( -9  >  ) 

28  PRINT  804.  I ,T,R .V.PRES.D.U.P.X.E 

804  FORMAT { 1M  . 1 4 , 1  E 1 2 . 4 . 1 F 1 2  •  7 . I E 1 4 . 5 •  1 F 1 4 . 5  •  A F 1 2 • 5 • 1 E 1 5 . 5 ) 

29  DT=DX/( <-A/T+V/R)*X) 

30  IF ( P ) 43 .31  .31 

31  T=T+DT 

32  R=R+V*DT 

33  U=U+DU 

3  4  P  =  P  f DP 

35  X  *  X  +  DX 

36  D  =  U  +  DD 

37  DCP=4.0*3.141593*P*D*CX*«4)*DX/IG-1.0) 

38  DTK  =  4. 0*3. 141593* (U«U) *  D*  (  X**4  )  *DX/2 . 0 

39  EP=EP+DEP 

40  EK~EK+DEK 

41  E=EP+LK 

42  CONTINUE 

43  GO  TO  1 

44  STOP 

4  5  END 
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NOMENCLATURE 


P 

P 

P1 

t 


r,R 


u 


V 

V 

V1 

7 

A 

a,p,6,e 

U,D,P,it,V,R 

P 

Eo 

kb 


s 

S 

C, 

M 

T 

r 


pressure 

density 

initial  density 

time 

radius 

radial  particle  velocity 
sound  speed 
specific  volume 
volume 

initial  volume 

adiabatic  exponent 
constant 

self -similar  exponents 
self-similar  functions 

y  ( 7-i)/(  7+1) 

total  energy 
kilobar 

arc  length,  constant 
entropy 

shock  wave  velocity 
C/c  Mach  No. 
temperature 
Griieneisen  factor 
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dynamic  process.  However,  actual  phenomena  show  that  the  irreversible 
zone  is  extremely  narrow. 

Outside  of  this  transition  zone  the  flow  obeys  the  laws  establish¬ 
ed  for  an  adiabatic  reversible  process.  Thus  the  empirical  fact  sug¬ 
gests  a  further  mathematical  idealization  i.e.,  an  infinitely  thin 
transition  zone.  We  find  that  the  original  differential  equations, 
valid  in  the  region  of  continuous  flow,  together  with  the  conditions 
expressing  the  conservation  laws  and  an  entropy  condition  across  the 
shock  front,  suffice  to  determine  the  flow  without  describing  in  de¬ 
tail  the  irreversible  process  across  the  shock  zone. 

There  is  considerable  difference  in  the  complexity  of  shock  waves 
in  gases  and  in  solids .  This  explains  the  greater  amount  of  litera¬ 
ture  on  the  former  compared  with  the  latter. 

The  difference  in  complexity  of  the  two  phenomena  is  attributable 
to  several  factors.  One  of  the  most  significant  being  the  lack  of  an 
equation  of  state  for  solids  which  is  valid  over  a  large  range  of 
pressures,  and  in  general,  a  lack  of  knowledge  of  the  thermodynamics 
of  solids.  (This  is  discussed  in  Section  4.9.)  Another  reason  is  the 
relative  difficulty  of  making  accurate  measurements  in  solids.  Other 
complications  are  (l)  a  shorter  time  sequence  in  solids,  (2)  ma¬ 
terial  imperfections,  (3)  non-homogeneity,  and  lack  of  transparency 
as  well  as  many  others. 

1.4.  Statement  of  Problem  and  Scope  of  Study 

The  problem  considered  in  this  research  is  that  of  the  dynamic 
behavior  of  a  ductile  thick  spherical  metal  shell  surrounding  a 
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(13) 

authors  and  employs  an  artificial  viscosity  as  a  mechanism  for 

avoiding  shock-front  discontinuities.  The  integration  process  con¬ 
sists  of  the  step-vise  solution  of  the  difference  equations  which  ap¬ 
proximate  the  differential  equations  of  motion  of  the  gas. 

Another  somewhat  different  approach  to  the  problem  of  a  strong 

(Ik) 

spherical  blast  wave  in  a  gas  is  that  used  by  P.  D.  Lax  '  Leix 

uses  the  conservation  form  of  the  hydrodynamic  equations  and  a  slight¬ 
ly  different  way  of  differencing  the  equations. 

One  of  the  first  investigators  of  the  cratering  problem  using 
numerical  techniques  was  Bjork.  In  predicting  the  crater  due  to  meter- 
oid  impact,  Bjork  assumed  the  impacted  earth  behaved  as  an  invis - 

cid,  compressible  fluid  and  postulated  an  equation  of  state  for  the 
material.  In  later  investigations  of  hypervelocity  impact  in  metals 
this  same  model  was  assumed  with  favorable  results 
Another  investigator  in  the  hypervelocity  impact  field  Is  J  M 

(17) 

Walsh  .  He  has  also  developed  a  numerical  program  for  the  inqoaet 
problem  based  on  the  "hydrodynamic"  model. 
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PUh2(|U’+U)  +  U(gU'+U)  +  (2P  +  gPD'/D  +  |P')  =  0  (3-12) 

(U-a)  gD'  +  6D  +  (gU'  +3tJ)D  =  0  (3-13) 

P'g(U-a)  +  P[2p-5(7-l)  +  2(U-a))-(D'/D)P|(7-l) 

(U-a)  =  0  (3.14) 

We  now  have  a  system  of  ordinary  differential  equations  for  the 
unknown  functions  U(g),  D(g),  P(g)  and  two  free  parameters  a  and 
5.  The  substitutions  of  Eqs.  (3-7),  which  may  appear  artificial,  are 
thus  justified. 

We  shall  now  reduce  the  number  of  variables  further.  Solving  Eq. 
(3.13)  for  gD'/i)  gives 

gD’/D  =  -  (&+|U'  +  3U)/(U-a)  (3-15) 

When  Eq  (3-15)  Is  put  into  the  remaining  two  equations  we  obtain  the 
pair. 


(u-a)(gu'+u)+  pu  +  p[2  (8+|u’+3u)/(u-a)]  +  gp"  =  0  (3.16) 

-(7-l)(&+gu'+3U)  -  2^(r-i)5-2  (u-a)-(u-a)gP'/p=o  (3-17) 

These  equations,  linear  in  gU'  and  gP'  may  be,  solved  simultaneous- 
ly,  giving 
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|U’  =  [-U(u-a)(U-l)  +  (&f2P»-3U7)P]/A 


(3-18) 


£P’  =  P  (7-l)u(U-l)  -  (U-a)[(37-l)u-2]  + 

((2g^2l3j6)  +  27)P  /A  (3.19) 

2 

where  A  =  (U-a)  -  7P. 

By  dividing  Eq.  (3»19)  by  Eq  (3  18) ,  the  following  ordinary 
differential  equation  is  obtained; 


dP  _  ill  P[H(U)  +  pq(u)] 
dU  |U’  =  R(U)  +  PfS(U)) 


(3  20) 


After  the  simplification, 

N(u)  =  7U(3a-l-2U)  +  (3-a)u-2a 
Q(U)  =  [2P  -(7-1)6]  /(U-a)  +27 

(3.21) 

R(U)  =  u(u-a)(l-u) 

S(U)  -  5+26+-  3U7 


This  is  the  basic  differential  equation  for  self -similar  motion.  After 
the  appropriate  solution  has  been  found  for  P  =  P(u),  the  function 
|  =  (u)  is  found  by  a  quadrature  of  Eq.  (3.18)  and  the  density  func¬ 
tion  D(|)  is  obtained  from  Eq.  (3.I5). 

As  we  shall  see,  these  self -similar  solutions  provide  a  suffic¬ 
iently  general  mathematical  description  of  an  expanding  cavity  reason¬ 
ably  consistent  with  the  given  conditions  of  iniation  of  the  process. 
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There  remains  the  problem  of  choosing  the  two  parameters  a  and  6. 

Let  us  assume  that  the  total  energy  given  up  by  the  explosive  ap¬ 
pears  as  kinetic  and  potential  energy  in  the  surrounding  medium.  This 
is  a  reasonable  one  for  the  cavity  expansion  process,  because  of  its 
short  duration,  provided  certain  secondary  effects  are  neglected.  With 
representing  the  shock  front  at  a  time  t,  the  total  energy  in 
the  fluid  shell  (potential  +  kinetic)  at  time  t  is  given  by 


=  /  (^T}  w2dr  +  / 


2  2 
pu  4jtr  dr 


(5-22) 


OL  Q£ 

where  rQ  =  |Qt  is  an  arbitrary  member  of  the  family  and  r1  =1^1 

is  the  location  of  the  shock  front. 

7 

The  first  term  arises  from  the  polytropic  relation  p=AP  which 
has  been  assumed  for  the  material  and  the  fact  that  the  work  done  by 
the  material  in  becoming  compressed  from  an  initial  volume  to  a 

a  volume  V  is  given  by 

V  p 

E(v)  =  f  pdV=P1V1  J'  Ap7(-^§)  =  -  715;+  const.  (5.23) 


where  we  have  used  the  relation  P2V2=P1V1‘  Usin6  the  substitutions 
in  (3.7)  the  energy  expression  becomes 
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P2(C-u2)  -  PXC  =  0  (329) 

P2u2(e-u2)  -  (p2-P;l)  =  0  (3.30) 

P2^Z2  +  7^1  ^~)(c_u2)  -  P2u2=0  (ei=°)  (3-31) 

It  is  useful  to  solve  Eq.  (3-29)  and  (3*30)  for  the  velocities.  Then 
u2  =  (p2-P1)/p1C  (3-32) 


If  p  is  in  kilobars,  p  in  gm/cc,  then  u  and  C  vill  come  out  In 
km/sec  as  expressed  by  the  following  formulas: 


U2  =  (p2  "  pi)/(10  P-lC) 


(3-34) 


It  is  useful  to  have  these  relations  for  a  polytropic  medium.  Then 
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u2  =  C  (1  -  n2) 

P2  =  P1  u2  C  =  P1  °2  ^ 

(3.36) 

p2^pl  =  (r+l)/(7-l)  =  ~2 

h 

°  -  M  -  A* 1,2  c 

The  last  quantity  is  conveniently  referred  to  as  the  "sound  speed"  in 
shock  wave  analysis.  These  relations  apply  Just  as  well  to  a  spheri¬ 
cal  or  curved  surface  as  to  a  plane,  since  the  effect  of  spherical  di¬ 
vergence  (the  2  u/r  term)  on  a  finite  or  sudden  jump  is  of  a  higher 
order  of  magnitude.  This  may  also  be  shown  geometrically  by  consider- 

-C£ 

ing  an  infinitesimal  surface  element  of  the  shock  front.  Since  |-j=rt 
(by  definition)  along  the  shock  front, 

-a  -a-l 

0  =  t  dr  -  cert  dt 

so  C  =  —•  =  art-1  =  cee  ta  .  Then  from  Eq  (4.7)  and  the  relations 
du  1 

(3.29)  through  (3>3l)  with 

C  -  u2  =  ce^t0'1  -  tP61U(61)  =  tp(a-u). 


we  obtain 


t6+P  |x  D(a-U)-  Pl  a  tp  -  0  (3.37) 

t&+2P  ^2  m  (a.u)  _  DP  _  o 


(3-38) 


22 

Relations  can  then  be  determined  for  D=D(g),  U=U(|),  P=P(g) 

using  equations  (3-15),  (3°l8)  and  (3-19). 

A  discussion  of  the  singular  points  of  the  differential  equation 
(3-20)  and  the  solution  curve  (3° ^5)  are  given  in  Chapter  IV. 
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A  numerical  integration  of  the  differential  equation  (3*20)  was 
performed.  (See  Appendix  B).  At  the  same  time,  particle  "paths"  were 
traced  out  and  the  -variation  of  the  physical  variables  along  these  par¬ 
ticle  paths  were  calculated.  (See  Appendix  B). 

Pressure  isobars  were  then  drawn  connecting  the  points  of  equal 
pressure  in  the  medium  (See  Fig.  5). 

Particular  attention  was  paid  the  particle  which  started  at 
r  =  rQ  =  1  698  as  this  marked  the  cavity  surface  and  was  of  great  in¬ 
terest  in  the  study.  (Figs.  5,6, 8,9)-  It  is  also  convenient  to  tabu¬ 
late  the  energy  integral  (3-24)  and  compare  it  at  some  selected  time 
to  the  energy  given  up  by  the  explosive 

In  the  theory  of  self -similar  motion,  an  assumption  of  constant 
energy  was  made  (See  Eq.  (3-24)  )  in  order  to  provide  the  condition 
(3  25)  for  determining  a  and  with  it,  all  the  other  exponents.  The 
energy  integral  (3.24)  is  extended  between  two  points,  one  of  which  is 

located  on  the  shock  front  £  =  E,  and  a  lower  value  £  =  £  •  The  in- 

al  o 

tegral  path,  such  as  BC  in  Fig.  6  may  be  arbitrarily  chosen,  so  long  as 
it  terminates  on  these  two  curves. 

The  energy  in  the  disturbed  part  of  the  solid  will,  however,  change 
with  time  because  its  lower  boundary,  the  cavity  surface,  is  not  one  of 
the  family  of  £-curves.  Thus,  the  energy  values  in  Table  (l)  repre¬ 
sent  an  integration  of  the  expression  in  (3-25)  taken  along  the  cavity 
surface  curve.  If  we  extend  this  integration  far  enough  (say  to  100 
microsec)  so  that  point  D  practically  coincides  with  point  C,  the 
values  in  the  table  become  asymptotically  constant,  and  we  have 


DIAGRAM  OF  SHOCK  REGION  FOR  ALUMINUM  SPHERE 
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TABLE  1 

Variation  of  Physical  Quantities  Along 
Solution  Curve  and  Cavity  Surface 


T 

U 

P 

1 

D 

E 

R 

sec  x 

ergs  x 

10-6 

gm/cc 

1012 

cm 

.8684 

.09302 

.02856 

451.28 

3.5182 

-0.0 

1.698 

•  9539 

.09312 

.02854 

438.45 

3.4282 

-.188 

1.713 

1.062 

.09324 

.02852 

424.25 

3.3282 

-  368 

1.730 

1.186 

-09337 

.02844 

410.10 

3.2282 

-.520 

1-748 

1.329 

.09352 

.02846 

396.02 

3.1282 

-.648 

1.767 

1.495 

.09369 

.02843 

381.99 

3.0282 

-755 

1.786 

1.688 

.09389 

.02840 

368.02 

2.9282 

-.844 

1.607 

1.939 

.09411 

.02836 

354.12 

2.8282 

-  918 

I.828 

2.181 

.09436 

.02832 

340.29 

2.7282 

-979 

1.851 

2.496 

.09466 

.02827 

326.52 

2.6282 

-1,029 

1.874 

2.872 

.09500 

.02821 

312.83 

2.5282 

-1,069 

I.900 

3.323 

•09539 

.02815 

299.21 

2.4282 

-1.101 

1.926 

3-868 

.09586 

.02807 

285-69 

2.3282 

-1,127 

1-954 

4-533 

.09641 

•02799 

272.24 

2.2281 

-1.148 

1.984 

5-350 

.09707 

.02789 

258.89 

2.1282 

-1.164 

2.016 

6.364 

•09787 

.02777 

245.63 

2.0282 

-1,176 

2.050 

7-636 

.09884 

.02764 

232.48 

1.9282 

-1.186 

2.087 

9-252 

.10003 

.02749 

219.43 

1.8282 

-1.193 

2.127 

11.33 

.10152 

.02731 

206.51 

1.7282 

-1.198 

2.171 

14.04 

.10341 

.  02710 

193.71 

1.6282 

-1,202 

2.219 

17.65 

.10585 

.02686 

181.05 

1.5282 

-1.205 

2.272 

22.52 

.10904 

.02669 

168.55 

1.4281 

-1.207 

2.332 

29.28 

•11331 

.02628 

156.21 

1.3282 

-1  209 

2.401 

58.91 

.11914 

.02593 

144.08 

1.2282 

-1.210 

2.481 

53.07 

. 12728 

•02555 

132.17 

1.1282 

-1. 211 

2.577 

74.76 

.13883 

.02513 

120.56 

1.0282 

-1.211 

2.696 

109.5 

.15536 

.02462 

109.38 

.9282 

-1.212 

2.849 

168.3 

.17866 

.02384 

98.85 

.8282 

-1 , 212 

3.058 

273.2 

.21002 

.02241 

89.36 

.7282 

-1.212 

3-355 

470.8 

.24867 

.01974 

81.30 

.6282 

-1.212 

3-795 

865.8 

.29068 

.01550 

74.96 

.5282 

-1.212 

4.465 

1717 

•32972 

.01012 

70.40 

.4282 

-1.212 

5.514 

3681 

•35941 

.00478 

67.43 

.3282 

-1.212 

6.968 

7404 

.37363 

.00098 

65.87 

.2282 

-1.212 

9-257 

8357 

.36027 

.00001 

65.67 

.1282 

-1.212 

9.687 

kk 


b.k.  Blast  Waves  in  Other  Metals 

The  study  of  blast  waves  in  other  metals  is  essentially  the  study 
of  the  part  7  plays  in  the  theory  of  self -similar  motion.  For  this 
reason  the  present  section  will  explore  in  some  detail  the  many  intri¬ 
cate  ways  that  7  enters  into  the  problem. 

Selected  for  study  were  3  common  metals,  copper  (7  =  85 ) >  iron 
(7  =  16),  and  lead  (7=5-9)  Figure  4  is  a  log -log  plot  of  the  shock 
wave  data  from  Ref.  19°  Straight  lines  were  passed  through  the  data 
and  the  slope  measured  to  obtain  the  values  of  7. 

In  the  method  of  obtaining  the  solution  of  the  differential  equa¬ 
tion  (3  20),  the  functions  U,  P,  D  and  the  particle  trajectories  as 
well  as  shock  wave  paths  are  the  same  as  In  Sec.  4.3-  All  that  is 
necessary  is  to  change  the  value  of  7  in  the  computer  program.  (See 
Appendix  B). 

For  purposes  of  comparison  the  distribution  of  velocity  and  pres¬ 
sure  were  calculated  on  the  cavity  surface  and  plotted  in  Figs.  8  and 
9-  An  r-t  diagram  for  the  four  metals  is  shown  in  Fig.  7- 

Of  particular  Interest  is  the  fact  that  the  7's  used  range  from 
5-9  for  lead  to  16  for  iron  This  value  used  for  iron  is  supposed  to 
represent  the  lower  pressure  range  in  which  iron  has  a  body  centered 
cubic  lattice.  It  can  be  shown  (Fig.  10)  that  for  7=7  the  shock 
front  falls  on  one  of  the  singular  points  of  the  differential  equation 
(3-20).  For  this  case  there  exists 
of  equation  (3. 20),  i.e.. 


a  particularly  simple  solution 
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(  ujo)  snipoa 


FIG.  7  -  CAVITY  EXPANSION  FOR  FOUR  METALS 
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have  the  same  values  at  the  shock.  This  allows  the  density  ratio 
p^/p^  to  vary  and  thus  accounts  for  one  of  the  non-similar  factors 
mentioned  earlier.  However,  even  though  we  have  satisfied  the  condi¬ 
tions  at  the  shock  front,  the  flow  field  behind  the  shock  is  still 
incorrect  For  the  cratering  problem  this  ie  shown  to  be  a 

fairly  good  approximation. 

We  recall  that  the  density  ratio  was  giver,  by  the  expression 

p2/pi  =  /-I 

This  can  be  rewritten  as 

Pg/ +  1 

7  =  p2/'p1  -  I  (4.25) 

Thus  if  the  shock  speed  were  such  as  to  cause  a  density  ratio  of  1,2, 
then  the  solution  is  assumed  to  be  the  self -similar  solution  for 
7  =  11. 

Note  that  as  Pg  ->  p^  the  value  of  y  becomes  infinite.  (As 
the  shock  wave  decays  to  a  stress  wave  y  -»  co  . ) 

If  the  shock  speed  Rg  is  considered  as  the  independent  variable, 
then  Eq.  (4.25)  becomes 

(2S  -  1)  R  +  c 

7  .  - - , - 2 -  .  (4.26) 

Rs  -  C 

The  relation  Rg  =  c  +  SUg  used  in  this  equation  has  been  shown 

(22) 

v  to  approximate  the  Hugoniots  of  several  metals  very  well. 
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Application  of  this  theory  to  the  problem  of  cavity  expansion  in 
metals  was  described  in  setae  detail  without  presenting  numerical  re¬ 
sults  . 

In  considering  the  application  of  this  theory  to  other  materials, 
the  only  difference  to  be  encountered  would  be  the  initial  value  of  y 
used. 
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detected  in  the  post -shot  measurements.  However ,  the  values  obtained 
from  the  recorded  data  presented  in  Table  3  do  reflect  this  non-uniform 
deformation.  The  values  in  Table  2  indicate  that  the  average  increase 
in  the  outer  diameter  of  the  small  spheres  is  about  2$  and  about  1$  for 
the  larger  spheres.  For  the  inner  cavity  the  post-shot  diameter  was 
found  to  be  about  75$  greater  than  the  pre-shot  diameter  for  the  small¬ 
er  spheres.  The  increase  for  the  larger  spheres  is,  however,  greater 
than  100$.  In  Table  3  we  see  the  values  of  the  wave  velocity  agree 
well  with  the  known  elastic  wave  velocity  in  the  aluminum.  It  is  also 
evident  that  the  free  surface  velocity  for  the  two  materials  is  apprec¬ 
iably  different. 

Figure  13  shows  a  typical  plot  of  data  showing  the  response  of  an 
18  cm  sphere  recorded  through  a  condenser  micrometer.  Notice  the  de¬ 
finite  tendency  toward  oscillation  of  the  outer  surface  Figure  ll*  is 
an  actual  trace  of  the  surface  response  of  the  25  cm  sphere  recorded  by 
the  Polaroid  camera. 

5-3-  Discussion  and  Conclusions 

One  of  the  most  significant  questions  that  we  hoped  the  experiments 
would  answer  was  how  much  the  cavity  size  was  influenced  by  the  loca¬ 
tion  of  the  outer  boundary  More  specifically,  how  large  the  sphere 
had  to  be  to  be  considered  an  infinite  medium  This  was  not  answered 
quantitatively,  but  it  was  found  that  the  location  of  the  boundary  was 
a  significant  factor.  We  notice  that  the  post-shot  cavity  was  consis¬ 
tently  larger  for  the  25  cm  sphere  when  compared  to  the  l8  cm  sphere. 
This  might  be  explained  through  a  coupling  effect  between  the  reflect- 
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flow.  For  this  reason  It  is  suggested  that  some  numerical  results  he 
obtained  using  the  theory  described  in  Section  4-5- 

In  order  to  obtain  an  exact  analytical  solution,  the  boundary  con¬ 
ditions  on  the  inner  cavity  surface  must  be  prescribed  for  all  values 
of  time.  This  will  require  a  detailed  analysis  of  the  effect  of  con¬ 
finement  on  the  detonation  pressures  of  the  explosive,  and  also  con¬ 
sideration  of  pulsation  of  the  exploded  gas. 

There  remains  the  problem  of  defining  the  equation  of  state  of  a 
solid  throughout  the  range  of  pressures  encountered  in  the  blast  wave 
problem.  This  problem  is  complicated  by  phase  changes  and  other  com¬ 
plicated  mechanical  behavior.  This  should  continue  to  be  an  area  of 
great  effort. 

There  are  several  questions  that  can  be  answered  by  experimental 
investigations  which  would  facillitate  the  formulation  of  a  satisfact¬ 
ory  theory. 

(1)  What  part  do  the  reflected  waves  play  in  the  final  cavity  size. 
This  could  be  answered  using  a  flash  x-ray  or  some  similar  device. 

(2)  What  is  the  actual  Initial  detonation  pressure  exerted  on  the 
metal. 

There  is  need  for  a  pressure  sensing  device  which  will  survive  the 
extreme  conditions  experienced  in  this  test. 

(3)  How  much  time  is  required  for  the  cavity  to  assume  its  final  con¬ 


figuration. 
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ing  and  the  density  was  decided  to  toe  a  better  parameter.  Using  the 
relation 


dD  =  -  d[(H^)u+_3U((U^)2  -  P)3  ..  . .  *» - 

y  (n+pq)2  p2  +  (r+ps)2 


a  program  was  developed  where  the  density  is  the  independent  variable 
and  a  physical  constraint  was  maintained  on  the  integration. 
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CHAPTER  I 
INTRODUCTION 

1.1  Reasons  for  Studying  Spherical  Shock  Waves  in  Solids 

The  study  of  spherical  shock  waves  in  solids  poses  an  interesting 
problem  both  from  the  standpoint  of  the  desirability  of  the  results  as 
well  as  the  inherent  theoretical  interest  of  the  problem. 

Information  gained  in  this  study  should  be  of  interest  to  investi¬ 
gators  in  many  fields .  The  dynamic  expansion  of  a  spherical  cavity  in 
a  metal  under  the  conditions  produced  by  a  high  explosive  is  closely 
related  to  several  important  problems.  The  phenomena  associated  with 
cratering  due  to  hypervelocity  impact,  auto-frettage,  underwater  ex¬ 
plosions,  underground  nuclear  testing,  nuclear  reactor  explosions,  and 
many  more  have  much  in  common  with  the  cavity  expansion  problem. 

The  study  of  shock  waves  in  solids  and  their  applications  to  the 
study  of  properties  of  solids  has  become  an  active  area  of  engineering 
and  physics . 

Propagation  of  these  intense  waves  through  solids  with  pressures 
up  to  9  megabars  has  yielded  fundamental  thermodynamic  data  which  could 
not  be  obtained  by  any  other  means.  Shocks  have  been  used  to  produce 
changes  in  crystal  structure  of  some  materials.  It  has  also  been  found 
that  shocks  can  change  electrical  insulators  into  conductors.  Shocks 
harden  metals,  create  and  alter  vacancies  and  dislocations  in  lattices 
and  shock  induced  impact  bonding  of  metals  is  now  commercial. 

It  is  felt  that  the  simple  geometry  and  relative  ease  of  perform¬ 
ing  experiments  warrants  a  rigorous  study  of  spherical  shock  waves  in 


solids . 
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scientific  and  engineering  people.  The  label  shock  wave  is  often 
attached  to  both  strong  and  weak  discontinuities. 

Weak  discontinuities  are  fronts  across  which  the  time  and/or  space 
derivatives  of  the  physical  variables  which  describe  the  motion  and 
state  of  the  medium  may  experience  a  finite  change.  These  will  be  re¬ 
ferred  to  as  "sound  waves".  These  waves  are  the  characteristics  of  the 
linear  differential  equations  that  describe  the  motion.  Across  a  sound 
wave  the  physical  variables  experience  only  infinitesimal  changes. 

Strong  discontinuities  are  fronts  across  which  the  physical  var¬ 
iables  which  describe  the  motion  and  the  state  of  a  medium  experience 
a  finite  change . 

There  are  two  types  of  discontinuity  surfaces,  contact  surfaces 
and  shock  fronts.  Contact  surfaces  are  surfaces  separating  two  parts 
of  the  medium  without  any  flow  across  the  surface.  Shock  fronts  are 
discontinuity  surfaces  across  which  mass  flows. 

For  these  waves  there  exists  no  sonic  transmission  of  the  discon¬ 
tinuities  of  the  physical  variables  since  the  governing  differential 
equations  are  nonlinear  We  find  that  the  discontinuities  of  the 
physical  variables  cannot  be  described  by  a  passage  to  a  liinit  from 
continuous  solutions  as  in  the  case  of  sound  waves  To  arrive  at  an 
adequate  theory  the  physical  fact  of  the  existence  of  the  shock  wave 
must  be  taken  into  account 

It  cannot  be  assumed  that  the  motion  is  isentropic,  since  the 
velocity  and  temperature  gradients  are  not  small.  Thus  we  are  con¬ 
fronted  with  having  to  describe  mathematically  an  irreversible  thermo- 
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Dli2  (|  U2  +  -^  P)  ^  (o-U)  -tBf3Pl15  DFU  =  0  (3  39) 


We  note  that  the  time  factor  cancels  in  Eq.  (3.38)  and  Eq.  (3-39)>  so 
that  these  equations  are  automatically  satisfied.  However,  to  secure 
independence  of  time  in  Eq  (3»37),  it  is  necessary  tp  make 


8-0  .  (3.40) 

With  this  condition  and  the  relations  (3.25),  the  assumed  form  for  the 
progressing  wave  solutions  reduce  to 

u  =  (r/t)  U(|)  P  =  (r/t)2  D( | )P( 5 ) 

P  =  D(|)  p/p  =  (r/t)2  P(|)  (3.41) 

with  |  =  rt”0 


This  solution  shows  that  on  the  shock  front,  where  |  is  constant,  the 
physical  quantities  such  as  velocity,  pressure,  density,  and  wave  vel¬ 
ocity  are  constant  on  the  rays  r/t  =  constant  This  also  dimensional - 

izes  the  functions  (4.7)  correctly,  for,  with  4  having  the  dimensions 

,  .a 

of  length  per  v^ime)  ,  and 

4  =  [I/T-0]  D(|)  =  [ML-5] 

U(4)  -  Cl]  HO  =  til 


then  u,  p,  P, 
ocity,  pressure 
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Another  means  of  approaching  the  problem  is  to  consider  a  flew 
field  which  depends  on  a  non-dimensional  parameter  T)  «  r/RQ(t) .  At 
each  instant  the  distributions  of  velocity,  pressure  and  density  are 
taken  to  be  the  same  when  viewed  on  a  scale  defined  by  the  shock  radius 
at  that  instant.  Then  if  we  assume  the  radius  of  the  shock  wave  fol¬ 
lows  a  law  Rs("k)  =  A  tn,  we  find  we  have  essentially  the  same  formu¬ 
lation  as  discussed  earlier.  All  that  remains  is  the  introduction  of 
the  dimensionless  velocity,  pressure  and  density  functions  which  are 
functions  of  tj. 

3-2.  Solution  of  the  Governing  Differential  Equation  of  Self -Similar 
Motion 

The  analytical  solution  of  the  governing  ordinary  differential 
equation  is  obtained  in  an  indirect  way  (see  Sedov  (Ref.  4)  ).  We  re¬ 
call  that  in  the  first  section  of  this  chapter  we  used  the  fact  that 
the  energy  integral  was  time  independent  to  obtain  a  relationship  be¬ 
tween  the  parameters  a  and  5 . 

It  can  be  shown  that  the  energy  assumption  implies 


!5[HJ 


(U-a)  (z|^  + 


constant* 


(3*44) 


This  must  be  true  for  all  stages  of  the  process,  in  particular 
on  the  shock  front.  Evaluation  of  the  constant  from  the  Ranklne- 
Hugoniot  conditions  shows  it  to  be  zero.  Hence,  we  have  effectively 
obtained  an  integral  of  Eq.  (3*20)  and  the  relationship  between  P  and 


U  becomes 


P 


(7-1)  u2  (u-a)/2  (a-7U) 


(3  45) 
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CHAPTER  IV 

APPLICATION  OP  THE  THEORY  OF  SELF -SIMILAR  MOTION  TO 
THE  PROBLEM  OF  EXPANSION  OF  A  SPHERICAL  CAVITY  IN 
A  METAL  DUE  TO  EXPLOSIVE  IMPACT 
4.1  Description  of  Model  and  Assumptions 

Detonation  of  a  high  explosive  confined  within  and  in  contact  with 
a  ductile  metal  generates  a  strong  shock  wave  which  raises  the  pressures 
abruptly  (i.e„,  in  negligible  time)  into  the  megabar  range.  These 
pressures  are  large  compared  to  the  material  strength  even  at  high 
strain  rates  and  one  is  led  to  assume  that  the  medium  behav<  3  as  an  in- 
viscid,  compressible  fluid- 

Thus  the  solution  of  the  problem  of  spherical  blast  waves  in  a 
metal  originating  at  a  spherical  cavity  reduces  to  solving  the  equa¬ 
tions  of  fluid  mechanics  (3-1)  -  (3-4)  (i.e.,  conservation  of  mass, 
momentum  and  energy,  together  with  the  equation  of  state  of  the  medium). 

If  we  next  assume  that  the  motion  is  self -similar,  then  all  of 
the  ideas  discussed  in  Chapter  III  are  at  our  disposal  The  remaining 
question  is  that  of  the  equation  of  state.  In  this  work  it  has  been 

7 

assumed  to  be  of  the  form  p  =  Ap  .  This  is  justified  to  some  extent 
in  the  next  section. 

4.2.  Equations  of  State  of  Metals 

The  equation,  of  state  plays  a  central  part  in  the  application  of 
the  theory  of  self -similar  motion  to  shock  waves  in  metals.  In  ad¬ 
dition  to  furnishing  a  relationship  between  pressure  and  density,  it 
furnishes  a  criterion  for  the  termination  of  the  regime  in  which  the 
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& 


for  a  QPT 

c2  =  ~  (Ap7  -  B)  =  7  Ap7'1 
For  later  use  we  will  define 

—  V 

p  =  p  +  B  =  Ap  . 


Consider  next  specific  internal  energy.  For  a  FT  medium; 


For  a  QPT  medium; 

e  =  J  -  (Ap7  -  B)  d(i) 

(  _A_  Bp~\ 

•  7-1  P  -1  ' 

=  (  +  B)V 

e  -  (E^jZSjV  <*-6) 


The  Hugoniot  function  is  defined  by 


P  + 


H(V,p)  =  e(V,p)  -  e(VllPl)  +  (V-VjK  g  -) 


(^•7 ) 


If  we  set  H(V,p)  =  0,  the  Hugoniot  relation  is  obtained.  This  rela¬ 
tion  characterizes  all  pairs  (V,p)  for  the  state  on  one  side  of  the 
shock  front  that  are  compatible  with  the  three  Si ock  front  conditions 
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also 


pa 4  M. E 

2  2  7 

p  p2  +  nl 


P2  +  ^  pl 
P1  +  P2 


(4.9) 


Thus  we  see  that  for  a  OPT  medium  the  ratio  —  is  not  a  constant  as 

P1 

in  the  case  of  a  FT  medium;  since  /  0.  Next  consider  using  the  QPT 
relation  for  self -similar  motion  In  order  to  obtain  the  solution 
curve  of  the  differential  equation  we  must  have  a  unique  starting  point 


in  the  P-U  plane.  Since  the  ratio  —  is  not  constant  there  exists  no 

P1 


such  point  and  we  are  forced  to  abandon  the  QPT  relation, 
it-. 3-  Solution  of  the  Blast  Wave  Problem  in  Aluminum 

For  a  shock  wave  moving  into  polytropic  medium  at  rest,  the  Ran- 
kine-Hugonlot  equations  become 


P2  (c  -  ug)  -  Px  C  =  0 

(4.10) 

P2U2  (C  '  U2^  '  ^p2  "  Pl^  = 

0 

(4  11) 

,1  2  1  P2w_  v 

p2(2  ug  +  JZx  p~^C  ~  u2' 

'  P2U2  =  ° 

(4.12) 

(ex  =  0) 

It  is  interesting  to  estimate  the  accuracy  of  the  assumption  of  a  strong 

shock  wave. 

To  do  so  the  above  equations  are  rewritten  in  another  form  using 

(3.32)  through  (3.36).  That  is 
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We  new  turn  our  attention  to  the  energy  integral  Eq.  (3- 25). 
R 


E  =  J  (^)+  \  u2)  P^Kr2dr 


=  4*  f  [(|)2f?i  +  \  (f)2  4-2dr.  (4.27) 


hut 


§.  =  R  t 
*1  s 


~  dR  .  R  .  _  _  R 

.  sr  -  <V  ■  <^(-|> 


also 


2  2  2a 

r  =  |  t 


Using  these  relations  and  substituting  where  appropriate  into  Eq. 
(4.25),  we  find 


E  =  R  2  Rq3  F(7)  , 


s  s 


(4.28) 


where 


<Sr  *  I”2)  £**  • 


Solving  (4.27)  for  R  we  find 

s 


R  =  ( - Vg) 

F(7)  Rs2 


1/3 


(4.29) 


Since  7  is  a  function  of  Rg  by  Eq.  (4.26)  this  expression  becanes 
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different  state  of  the  material,  about  which  information  is  almost  com¬ 
pletely  lacking.  Several  mechanisms  have  been  suggested  for  terminat¬ 
ing  the  cavity  expansion: 

(1)  an  energy-level  criterion 

(2)  a  temperature  criterion 

(3)  a  yield -point  criterion 

Criteria  such  as  (2)  or  (3)  are  tempting  because  they  tend  to  pro¬ 
vide  fairly  definite  marks  as  to  when  the  material  "freezes",  either 
when  a  given  temperature,  or  a  given  pressure  is  reached.  However,  our 
knowledge  of  materials  under  these  conditions  is  still  too  incomplete 
to  solve  this  problem.  The  total  energy  of  the  moving  material  appears 
to  stop  increasing  after  the  expansion  phase  has  begun,  so  there  is  no 
change  in  energy.  Furthermore,  any  quantitative  use  of  energy  balances 
would  require  careful  accounting  of  all  the  energy  losses  as  well. 

It  was  shown  that  in  the  self -similar  approach  the  only  difference 
in  the  approach  used  for  different  materials  was  the  material  constant 
7.  Also  there  is  considerable  difference  in  the  behavior  of  materials 
whose  value  of  7  is  less  than  7  compared  to  those  having  a  value  of 
7  greater  than  7>  due  to  the  completely  different  behavior  of  the 
solution  curve.  (See  Fig.  10)  One  must  be  reminded  that  these  compari¬ 
sons  are  valid  only  under  the  assumption  that  the  equation  of  state 
used  in  this  analysis  is  valid  in  the  stated  range  of  pressures. 

In  the  previous  section  a  procedure  was  discussed  for  extending 
the  regime  in  which  a  self -similar  type  solution  is  valid. 
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Specimen  Failed 
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The  permanent  deformation  of  the  inner  and  outer  surface  of  the 
sphere  as  veil  as  the  material  damage  vere  of  primary  concern  in  the 
experiment.  Hcrvever,  there  vere  other  quantities  that  vere  of  great 
interest  such  as  free  surface  velocity  time  curves  3  time  of  arrival  of 
first  disturbance  and  maximum  deflection  and  shock  velocity  time  curve 
througi  the  material. 

Two  condenser  type  micrometers  were  used  to  obtain  the  free  sur¬ 
face  velocity  and  maximum  radial  deflection.  These  micrometers  consist 
of  a  condenser  arrangement  in  which  the  surface  of  the  sphere  consti¬ 
tutes  the  ground  side  and  the  other  plate  is  placed  a  few  tenths  of  a 
centimeter  away.  (See  Fig.  12).  The  output  of  these  micrometers  is 
recorded  on  cathode  ray  oscilloscopes  which  are  triggered  by  a  simple 
ionization  probe  inserted  through  the  plug  into  the  explosive  detonator 
cavity. 

Quartz  disc  crystals  are  affixed  to  aluminum  rods  and  threaded  in¬ 
to  the  sphere  to  varying  depths  to  measure  the  time  of  arrival  of  the 
disturbance  through  the  aluminum  sphere.  The  output  from  the  crystals 
is  recorded  in  the  same  way  as  the  condenser  micrometers . 

After  the  apparatus  is  assembled  the  condenser  micrometers  are 
calibrated  remotely.  The  explosive  is  then  initiated  and  the  scope 
traces  recorded  with  still  cameras  using  Polaroid  film. 

5.2.  Test  Results 

Table  2  contains  a  tabulation  of  pre-shot  and  post -shot  physical 
measurements  of  each  specimen.  The  inner  cavity  and  outer  surface 
measurements  are  averages  and  do  not  reflect  the  asymmetric  deformation 
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FI6.  12  EXPERIMENTAL  ARRANGMENT 
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the  metal  is  fairly  well  established,  but  terminates  under  another, 
about  which  there  is  little  information. 

In  applying  the  theory  of  self-similar  motion  to  cavity  expansion 
in  other  metals,  it  is  seen  that  the  only  difference  is  the  material 
constant  7  The  behavior  of  the  solution  curves  for  variable  7  is 
shown  in  Fig.  l6  For  7=7  the  solution  is  particularly  simple  (see 
Eqs  ,  (4.24)  ). 

In  considering  the  application  of  the  theory,  which  takes  into  ac¬ 
count  the  variable  density  ratio  across  the  shock  front  (see  Sec.  4-5)> 
we  see  that  a  numerical  approach  can  be  successfully  used  to  extend  the 
self -similar  theory. 

In  Chapter  V  the  results  of  new  experiments  were  presented.  These 
studies  were  made  for  the  purpose  of  verifying  the  foregoing  theory- 
Aluminum  spheres  of  18  cm  and  25  cm  diameter  were  subjected  to  an  in¬ 
ternal  blast  from  a  32  gn  spherical  charge  of  Pentolite.  The  result¬ 
ing  surface  motion  and  permanent  deformation  were  presented  and  the  ef¬ 
fect  of  the  size  of  the  sphere  established. 

The  inner  cavity  of  the  smaller  spheres  was  deformed  from  an  in¬ 
itial  diameter  of  J.kO  cm  to  6.0  centimeters.  The  cavity  of  the  larg¬ 
er  spheres,  however.  Increased  from  a  diameter  of  J.k  to  7-b  cm. 

Other  pertinent  data  is  given  in  Tables  2  and  5- 

The  actual  cause  of  the  difference  in  cavity  size  is  not  known. 
However,  the  possibility  of  coupling  of  the  waves  due  to  reflection 

from  the  boundaries  is  a  possible  explanation  Another  possible  cause 

(24) 

is  that  presented  by  Hill  '  in  which  the  motion  is  assumed  quasi- 
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dP  rs(P)  +  m  (£&+  2r(\J  -  a))] 
dU  U(l-U)  +  (2P  +  3P7te 


(A.6) 


We  now  seek  a  solution  curve  which  approaches  the  point  U=a  with 
a  definite  limiting  slope  (i.e.  be  locally  a  straight -line) .  Letting 
U  -»  a  in  Eq.  (A .6),  placing  the  resulting  egression  in  (A.  A),  and  di¬ 
viding  by  (U-Ct),  we  obtain  the  following  equation  for  m: 


[N(U)  +  m  (2P  +  2 7  (U  -  a))] 

m  =  m  U(l-U)  +  \%  l  3Uy;m - ^  (A*7) 

Equation  (A- 7)  has  the  following  three  solutions: 

1)  m  =  go 

2)  m  *»  0 

3)  m  -  (a-l)/5 

Thus  there  are  three  solution  curves  which  pass  through  the  point.  The 
first  two  are  the  straight  lines  U  =  cr,  P  =  0  respectively,  which 
are  singular  to  the  differential  equation.  The  remaining  slope  is 
m  =  -  0.2  for  Of  =  0.4. 

By  proceeding  along  the  straight  line  of  slope  m  a  short  dis¬ 
tance  toward  S4  and  using  this  point  as  the  starting  point  for  the 
numerical  calculation,  a  very  close  approximation  to  the  actual  curve 


is  obtained. 


